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Abstract: We argue that in a nonlinear gravity theory (the Lagrangian being an arbitrary function of 
the curvature scalar R), which according to well-known results is dynamically equivalent to a self-gravitating 
scalar field in General Relativity, the true physical variables are exactly those which describe the equivalent 
general-relativistic model (these variables are known as Einstein frame). Whenever such variables cannot 
be defined, there are strong indications that the original theory is unphysical, in the sense that Minkowski 
space is unstable due to existence of negative-energy solutions close to it. To this aim we first clarify 
the global net of relationships between the nonlinear gravity theories, scalar-tensor theories and General 
Relativity, showing that in a sense these are "canonically conjugated" to each other. We stress that the 
isomorphisms are in most cases local; in the regions where these are defined, we explicitly show how to map, 
in the presence of matter, the Jordan frame to the Einstein one and backwards. We study energetics for 
asymptotically flat solutions for those Lagrangians which admit conformal rescaling to Einstein frame in the 
vicinity of flat space. This is based on the second-order dynamics obtained, without changing the metric, 
by the use of a Helmholtz Lagrangian. We prove for a large class of these Lagrangians that the ADM energy 
is positive for solutions close to flat space, and this is determined by the lowest-order terms, R + aR^ , in the 
Lagrangian. The proof of this Positive Energy Theorem relies on the existence of the Einstein frame, since 
in the (Helmholtz-) Jordan frame the Dominant Energy Condition docs not hold and the field variables are 
unrelated to the total energy of the system. This is why we regard the Jordan frame as unphysical, while 
the Einstein frame is physical and reveals the physical contents of the theory. The latter should hence be 
viewed as physically equivalent to a self-interacting general-relativistic scalar field. 
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1. Introduction 

Metric theories of gravitation which are based on a Lagrangian density depending in a nonlinear way 
on the scalar curvature - such theories, somehow improperly, are usually called "nonlinear gravity" (NLG) 
models - share a general property, which has been extensively described in several works: acting on the 
metric by a suitable conformal transformation, the field equations can be recast into Einstein ones for the 
rcscalcd metric, interacting with a scalar field. It is therefore claimed that any NLG theory is equivalent to 
General Relativity (with the scalar field). 

A similar phenomenon occurs in Jordan-Brans-Dicke theory and its generalizations, the scalar-tensor 
gravity theories (STG). The original pair of variables (metric + scalar field) can be replaced by a new pair in 
which the metric has been conformally rcscalcd, and in the new variables the field equations become those of 
General Relativity; the scalar field, which is not affected by the transformation, turns out to be minimally 
coupled to the rescaled metric. The original set of variables is commonly called Jordan conformal frame, 
while the transformed set, whose dynamics is described by Einstein equations, is called EinsteAn conformal 
frame. A problem thus arises, whether the tensor representing the physical metric structure of space-time 
is the one belonging to the Jordan frame or to the Einstein frame. 

This problem can be traced back to Wolfgang Pauli in the early fifties (quoted in [1]). In a system 
consisting of metric gravity and a scalar field there is an ambiguity: the metric tensor can be conformally 
rescaled by an arbitrary (positive) function of the scalar. Thus, besides the original (Jordan) and Einstein 
frames, there exists an infinite number of conformally-related frames, each consisting of a pair (i^((/3)(7^^, (p) 
with a different F. One asks: are the metrics in the frames {Qhvt^) and {F{(p)gfj_^, f) also physically equiv- 
alent? The same question would clearly arise also for a more general change of variables g'^j^^ = g'^^{gai3, f), 

NLG and STG theories are in fact deeply connected, as is described in this paper, and can actually be 
viewed as different versions of the same model. We now define the notion of Jordan frame and Einstein frame 
for NLG theories; the main purpose of this paper is to analyze the problem of determining which "frame" ^ 
is the physical one. By "physical frame" we mean a set of field variables which are (at least in principle) 
measurable and satisfy all general requirements of classical field theory, e.g. give rise to positive- definite 
energy density (we are aware of the fact that the term "physical" is frequently abused in the literature). 

We assume, for simplicity, that the space-time is four- dimensional, although all the calculations can 

be actually carried out in higher dimension without significant modifications. The signature is ( — h ++) 
and we set ?i = c = SttG = I. Let us consider the Lagrangian of a vacuum NLG theory, 

L = f{R)^, (1.1) 

which generates the fourth-order equations 

f'{R)R^. - ^f{R)9^. - V^V./'(i?) + g^. □ f'{R) = 0, (1.2) 

where f'{R) = ^ and □ = 5'"^V^Vi/. In the vacuum case, the Jordan frame includes only the metric 
tensor (/^^. According to a well-known procedure [2-4], we introduce a pair of new variables {gnv,p), related 
to g^v (and to its derivatives) by 

P = f'{R), g^iu = Pg^iu- (1-3) 

(this transformation was rediscovered several times and generalized to the case of Lagrangians depending 
also on □ ^ii, [5, 6]) The scalar p is dimensionless, and to ensure the regularity of the conformal rescaling it 
is usually assunuxl that /) > (wc will retake this point below). Then the two metrics g^i, and ^j^i/ have the 

^ We use here the word "frame" to denote a choice of dynamical variables, rather than a choice of a 
reference frame in space-time; this might seem misleading, however this terminology is commonly used in 
the previous literature on the subject. Sometimes the term "gauge" (also abused in this context) is used 
instead of "frame" [10]. 
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same signature. Let r{p) be a solution of the equation ,f'[r{p)] —p = (it may not be miique, and we will 
consider this problem later); the field equations for the new variables become the following (second-order) 
ones: 

Gi^u =P'^ - \9ni^9'^^P,cxP,i3) + ]^{f[r{p)] - P • r-(p)),9^^| , (1.4a) 

whereby G is the Einstein tensor of the rescaled metric ^^j/, and 

Up-p-' ^r''P,,.P,u + \ (2/[r(p)] -pr(p))| = 0. (1.46) 
These equations can be derived from the Lagrangian [3] 

L = R^g - ^p-^r^P^^P,. +P-V(p) -p-V[r(p)]| V^- (1.5) 
It is customary to redefine the scalar field by setting p = eV?"^; the Lagrangian then becomes 

l=[r- r^M. - ^/^■, (1.6) 

one sees that is minimally coupled to and is an ordinary massive self-interacting scalar field with a 
potential 

V{ct>) = \e-^f r\p{cl>)] - ^e-'Vi'Pf(r\p{cP)]) , (1.7) 

which is determined by the original Lagrangian (1.1). The variables {Q/j.^, 4>) provide the Einstein frame for 
the NLG theory. 

In the Jordan frame, gravity is entirely described by the metric tensor g^y. In the Einstein frame, the 
scalar field </> acts as a source for the metric tensor g^^ and formally plays the role of an external "matterfield" ; 
however, the original theory did not include any matter, hence we are led to regard the scalar field occurring in 
the Einstein frame (which corresponds to the additional degree of freedom due to the higher order of the field 
equations in Jordan frame) as a "non-metric" aspect of the gravitational interaction. From this viewpoint, 
the NLG theory, although mathematically equivalent to General Relativity, is physically different because in 
the Einstein frame, where the equations coincide with those of GR, gravity is no longer represented by the 
metric tensor alone. However, another viewpoint is possible: one could assume that the fourth-order picture 
in Jordan frame represents an "already unified" model including a non-gravitational degree of freedom (a 
minimally coupled nonlinear scalar field), and that the gravitational interaction is described only by the 
rescaled metric in the Einstein frame (see Appendix D). 

This clearly indicates that the problem of the physical nature of the fields occurring in both frames 
should be addressed prior to any other consideration on the physical significance of the equivalence between 
NLG and GR. From the outline presented above, one might be led to the following conclusion: if one starts 
from a vacuum, NLG theory, there is no way to decide a priori which frame should be taken as the "physical" 
one: the choice of the physical metric is an additional datum, which affects the physical interpretation of 
the model but is essentially independent of the formal structure of the theory. On the other hand, if we 
formulate a gravitational theory including from the very beginning matter fields, the ambiguity is broken by 
the coupling of the metric tensor with such matter fields. In fact, the two metric tensors, g^i, and g^^, interact 
in a different way with each external field (we will discuss this point in greater detail below), therefore one 
should be able to single out the physical metric by requiring that matter fields be minimally coupled with 
it, and that neutral test particles fall along its geodesic lines. 

This viewpoint, more or less explicitly formulated, seems to be shared by many authors. However, a 
criterion based on the coupling of matter with gravity would be effective only if the total Lagrangian were 
determined on independent grounds, e.g. by string theory. As a matter of fact, this is not so: the theories 
which we consider in this paper are commonly viewed as fundamental ones and constructed by adding usual 
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interaction Lagrangians to a purely gravitational one. The point, that some authors seem to overlook, 
is that adding minimal-coupling terms to the Lagrangian (1.1) already entails that the Jordan frame is 
assumed to be the physical one. In this situation, claiming that the Einstein frame is unphysical, because 
the coupling between matter and the rescaled metric is not the usual one, is a petitio principii. It would be 
equally reasonable to add standard interaction terms, minimally coupled to g^i,, to the Lagrangian (1.6): the 
resulting coupling with gi^jy would then turn out to be unphysical. In other terms, the frequently repeated 
argument [7] in favour of the Jordan frame, namely that physical measurements arc made in this frame 
("atomic frame") since atomic masses are physical constants there, can be equally well used in favour of the 
Einstein frame, because the argument is a direct consequence of the arbitrary (at this level of reasoning) 
choice of the full action for a gravitational theory. Thus, the ambiguity is still present. 

Indeed, this ambiguity is faithfully represented in the literature. The authors dealing with nonlinear 
or scalar-tensor theories of gravity in the context of cosmology or of high-energy physics can be broadly 
divided into four groups. According to the authors of the first group, the Jordan frame is the physical one 
and the Einstein frame merely serves as a practical computational tool^[7, 8, 9]. (Barrow and Maeda [7] 
remark that the Einstein frame is computationally advantageous only in a vacuum theory. If matter, e.g. a 
perfect fluid, is included, the conformal transformation usually does not lead to simpler equations since 
in the Einstein frame the scalar field is coupled to the fluid). The authors in the second group regard the 
Einstein frame as being physical, either because of its resemblance to General Relativity"^ [10, 11], or since the 
standard formalism for quantizing the scalar field fluctuations in the linear approximation does not apply in 
the Jordan frame (or at least is suspect there) [12, 13], or because the masslcss spin-two graviton in Jordan- 
Brans-Dicke theory is described by the Einstein-frame metric [14], or finally because this is implied by 
dimensional reduction of a higher-dimensional action [15]. The third group consists of the authors claiming 
that the two frames are physically equivalent, at least at the classical level, since conformal transformations 
do not change the mass ratios of elementary particles and therefore does not alter physics'' ("physics cannot 
distinguish between conformal frames") [16]. The last, rather inhomogeneous group, involves authors who 
either use both conformal frames without addressing the problem of which of them (if any) is physical, or 
work exclusively in the original Jordan frame without making any reference to the existence of the rescaled 
metric [17]. 

In our opinion, the strongest argument which exists in the literature in favour of one of the conformal 
frames (if the universe is exactly four-dimensional) is the one based on quantization of field fluctuations. "It 
appears as if the quantization and conformal transformation are two mutually noncommutable procedures" 
[13] . Here we show that classical relativistic field theory provides another argument which points to the same 
conformal frame. 

The main arguments presented in this paper can be outlined in the following statements: 

(i) In contrast with claims raised in the previous literature, nonlinear and scalar-tensor theories of gravity 
can be equally well formulated, as far as formal consistency is considered, in terms of either one of the 
frames; from this viewpoint, either frame can in principle be assumed to be the physical one; 

(ii) However, the physics described by nonlinear and scalar tensor theories of gravity is not conformally 
invariant and using different conformal frames one finds inequivalent effects. Therefore, choosing a 
correct (i.e. physical) frame is an indispensable part of theoretical investigation [18, 19]; 

^ However, Kalara et al. [7], while investigating the power-law inflation in the Einstein frame, interpret 
the solutions and fit their parameters as if this frame were the physical one. 

^ Gibbons and Maeda [10] admit that there might be an argument (or merely a someone's prejudice) 
stating that the "physically correct" frame is different from the Einstein one. 

^ Garay and Garcia Bellido in [16] introduce a concept of "physical frame" which is different from ours. 
According to them the "physical" frame is one "in which observable particles have constant masses, since in 
this frame particles follow geodesies of the metric" . Then, by the assumption on the form of the full action, 
the Jordan frame coincides with the "physical" one. We notice at this point that in General Relativity 
particles (with constant masses) usually do not move along geodesies: consider for instance a perfect fluid 
with pressure. Geodesic motion is rather exceptional, it occurs in the case of pressureless dust and for 
non-interacting test particles. 
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(iii) The physical metric should be singled out already in the vacuum theory; the coupling of a given metric 
to matter fields is in fact determined by the physical significance ascribed to it, i.e. by its relation to 
the physical metric. 

(iv) In particular, if the physical metric is assumed to be the rescaled (Einstein-frame) one, the original non- 
linear Lagrangian should include interaction with ordinary matter in such a way that the corresponding 
coupling with the Einstein-frame metric turns out to be minimal. This problem has never been ad- 
dressed in the previous literature. We provide here a systematic and unique way to obtain Lagrangians 

which reduce to any given vacuum nonlinear Lagrangian in the absence of matter, and reproduce upon 
a suitable rescaling any (with some restrictions) minimally coupled matter Lagrangian; 

(v) Since consistency arguments do not allow one to establish which set of variables is physical, we investigate 
a different criterion. This criterion is purely classical and is provided by the positivity of energy for 

small fluctuations of every dynamical variable around the ground state solution; 

(vi) The notion of gravitational energy in NLG theories is inapplicable to this aim, due to the lack of a 
positivity theorem for higher order gravitational theories. To circumvent this problem, one should 
compare two second-order versions of the same theory, namely the second-order dynamics obtained in 
the Einstein frame and the "Hamiltonian" formulation of the theory constructed - without rescaling the 
metric - through a Legendre transformation; 

(vii) We show that the Positive Energy Theorem for Nonlinear Gravity, proven by Strominger [20] for a 
quadratic case, holds for a larger class of Lagrangians for which the Einstein frame can be defined 
around flat space. Actually, whether the ADM energy is positive or not depends on the potential (1.7). 
The existence of the Einstein frame is in any case essential for assessing classical stability of Minkowski 
space and positivity of energy for nearby solutions. In the Jordan frame, the Dominant Energy Condition 
never holds. For these reasons, the Einstein frame is the most natural candidate for the role of physical 
frame in NLG theories. 

The paper is organized as follows. In Section 2, we discuss in full detail the statements (i), (ii) and (iii) 
above, using examples taken from the literature. Some material is alrady known, nevetheless the amount 
of confusion existing in the recent literature on the subject justifies a detailed exposition. We provide 
there a global and complete picture of relationships between NGL and STG theories and general-relativistic 
models of a scalar field, while previously only separate connections were known. Section 3 is the heart 
of the paper. We discuss there items (v), (vi) and (vii). We investigate there a large class of nonlinear 
Lagrangians, for which flat space (in the Jordan frame) is a solution and the Einstein frame exists for it. 
While the ADM energy can be defined for any asymptotically flat solution, one is unable to establish its 
sign in the Jordan frame. In the Einstein frame, the relation betwen the interior of the system and its 
total energy takes on the standard form known from General Relativity. This fact convinces us that the 
Jordan frame is unphysical. Section 4 contains conclusions. All technical parts of the paper are moved to 
appendices and the main body of the paper can be read without consulting them In particular, Appendix 
C provides theoretical background for some results applied in Sections 2 and 3. The "inverse problem of 
nonlinear gravity" (finding the purely metric nonlinear Lagrangian which generates a prescribed potential 
for the scalar field) is presented in Appendix D, and the difficulties arising from possible local failures of the 
conditions ensuring the existence of the Einstein frame (a problem sometimes raised in the literature, but 
still lacking a systematic investigation) are dealt with in Appendix E. 
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2. Interaction of gravity with matter and conservation laws 

Some authors have tried to solve the problem of determining the physical metric by showing that only 
one possible choice allows one to obtain a divcrgcncclcss energy momentum tensor for any sclf-gravitating 
matter. Brans used this argument to claim that only the Jordan frame is physical [18], while in a recent 
paper [21] Cotsakis was led by a more detailed investigation of this point to the conclusion that the Eistein 
frame is the physical one. Unfortunately, an error invalidates Cotsakis' proof ^. Here wc show that studying 
the conservation laws for matter does not allow one to find out which frame is physical: the equations of 
motion for matter and gravity form consistent and closed systems for both Jordan and Einstein frames. 

Our basic assumption is that the gravitational interaction of matter should bedescribed by minimal 
coupling with the physical metric tensorfield. In other words, the physical metric of space-time should be 
identified prior to the construction of the Lagrangian for a gravitating system. Different identifications of the 
physical frame will give rise to (mathematically and physically) inequivalent Lagrangians. This assumption is 
not commonly accepted; e.g., the authors of [13] argue that the Einstein frame is physical, but nevertheless 
they assume that matter minimally couples to the metric in the (unphysical) Jordan frame. A similar 
viewpoint is taken by Alonso et al. in [11]. 

The assumption is based on the postulate that the groat advantage of Einstein's General Relativity, 
i.e. the universal validity of the minimal coupling principle, should be retained in NLG and scalar-tensor 
theories of gravity. The principle is partially abrogated in the so-called "extended Jordan-Brans- Dicke the- 
ory" , where the ordinary ( "visible" ) and dark ( "invisible" ) matter minimally couple to different (conformally 
related) metric fields (Damour et al. in [7]). Therefore, although observationally viable, we do not take this 
theory into account in this paper, on theoretical grounds. 

In the first part of this Section, we discuss the implication of our assumption in the context of NLG 
theory, considering separately the two cases: either the Jordan-frame or the Einstein-frame metric is regarded 
as physical. We show that the two cases describe different physical models of gravitational interaction of 
matter fields, and we show that each of the two models can be consistently formulated in both frames. 

The scalar-tensor gravity (STG) theories share with NLG theories the feature that they can be refor- 
mulated as a general- relativistic model for a self-gravitating scalar field^. As a matter of fact, any NLG 
theory can be also in terms of a scalar-tensor theory with nontrivial cosmological function (without changing 
the metric). This is explained in detail in Appendix C, and is extensively used in Section 3. Few authors 
seem to be aware of this relation, while many authors establish a connection between NLG and STG theories 
via the equivalence to General Relativity. Hence, NLG and STG theories appear strictly intertwined in the 
literature, sometimes causing confusion. A physical difference between STG and NLG theories is that in 
the former it is a priori postulated that matter couples minimally to the Jordan-frame metric. The dy- 
namical equivalence clearly shows that, despite appearances, there is no deeper difference between the two 
classes of gravity theories. Nevertheless, for clarity reasons, we also discuss here interactions with matter 
and conservation laws for STG theories. 

As it is shown in Appendix D, the "inverse problem of nonlinear gravity", i.e. finding a NLG Lagrangian 
which is equivalent to General Relativity plus a scalar field with any self- -interaction potential U {(f)) , has 
a solution. An exception is provided by the massless linear scalar field: if Uicj)) = and the cosmological 
constant vanishes, the equivalent nonlinear Lagrangian does not exist and the scalar field cannot be absorbed 
into the Jordan-frame metric. However, the massless scalar field is conformally equivalent to both a STG 
theory and to the conformally-invariant scalar field model. 

First we investigate the consequences of our assumption for the case of NLG theory. 



^ The crucial step of the proof consists in taking the divergence of the fourth-order equation (eq. (2.2) 
below): Cotsakis claims that it does not vanish in general. Actually, it is a generalized Bianchi identity (see 
Appendix A). 

^ It is worth noticing that we should discriminate between the dynamical equivalence of Lagrangians and 
equality of the action integrals, see Appendices C and E. This difference is significant in quantum theory; in 
classical field theory it is irrelevant as physical meaning is given only to solutions of the field equations. 
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2.1 NLG theory, case I: the original metric g^^, is physical. 

According to the assumption, the Lagrangian for gravity and matterfields (collectively denoted by reads'^ 

L = [f{R) + 2£^at(*, 9)]^/^■ (2.1) 
The gravitational field equations are then 

^[\/^/(^i)] = /'(i«)ii^.-^/(i?)5^.-V^V./'(ii)+5^.n/'(i?)=T^.(*,5) , (2.2) 

where, as usual, T^;^ = — ■y=-3^(\/— 5^mat)- It follows from Noether theorem (which in this case is 

equivalent to a generalized Bianchi identity) that Vj^T^'' = (sec Appendix A). Using the general procedure 
described in Appendix C and upon conformal rescaling (1.3) of the metric one gets the Einstein- frame 
Lagrangian for the system (metric + scalar field + matter) [3] 

L{g, cf>, f ) = [R{g) - g"''.^,^.^,. - 2V{<P) + 2e-2\^^Cat(f , e-^'^'g)] (2.3) 

which is dynamically equivalent to the Lagrangian (2.1). Here, as before, (p = y^lnp is a scalar field having 
the canonical kinetic term and self-interacting via the potential V{4>) given in (1.7). The field equations 
for the metric and the scalar in the Einstein frame can be obtained either by transforming equations (2.2) 
or directly from the Lagrangian (2.3). In the latter case, caution is needed whilst taking variations of the 
matter part, to avoid ambiguities in the definition of the stress tensor for the matter field '5. In fact, defining 
separate stress tensors for (j) and for "if makes sense only when the matter Lagrangian can be split in two 
parts, each one depending only on the metric and on one of the fields. Therefore, it is advisable to formulate 
the gravitational equations in both frames in terms of the tensor T^^ already defined in terms of the physical 
metric*. To this aim, one should first take the variation of the matter term with respect to * and to g'^'^ 
and then use Sg^"^ = p5g^" + g^'^Sp. As a result, one finds the field equations for the metric in the form 

G^, = t^4<f>, 9) + e-^%^{^, e-^f'g), (2.4) 



where ^/S = y 5 and 

tfiiy = 4>,^4i,v - -j^9fiug°'^^<i>.a<t>.p - y{4>)9pii' (2-5) 

plays the role of the effective stress-energy tensor for the scalar 0, and the equation of motion for (p 

y^^_dV 1 

with 

dV 1 



d(j) V6p 
The Bianchi identity now implies 



□ </> = ^ + 4e-VI^T^,(vI/,e-^/i^g) , (2.6) 

(p=p(</,)=eVT*) . (2.7) 



V6 

2 

-fHv)] - r{p) 



[V(0, 9) + e-\^^T^.(*, e-^*5)] = 0, (2.8) 

and the two stress tensors arc not separately conserved, since the scalar field interacts with matter, as is 

explicitly seen in Lagrangian(2.3). 

In this picture the scalar field influences the motion of any gravitating matter, except for the particular 
case in which the interaction lagrangian ly^at is conformally invariant (see Appendix B). Prom the physical 
viewpoint this theory is not equivalent to General Relativity: in fact, gravity is completely represented by a 
metric tensor only in the Jordan frame (where it obeys fourth-order equations), while in the Einstein frame, 
where equations of motion for g^,^ are formally Einstein ones, there is a non-geometric gravitational degree 
of freedom, represented by the scalar (/>, which is universally coupled to matter. From the viewpoint of 
conservation laws, however, no inconsistency arises from the assumption that the original metric g^i, be the 
physical one. 



^ The coefficient in front of ^mat is chosen on the assumption that the linear term in the Taylor expansion 
of f{R) has coefficient 1 (see Appendix A). 

* We notice that the definition provides the covariant components of the stress tensor; if contravariant 
components T^^"^ are defined by raising the indices with the physical metric g'^'^ , then the contravariant form 
of equation (2.4) reads G^"' = t^'^icj), g) + e'^'^'t'Ti^'' , e-'^'^g). 
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2.2 NGL theory, case II: the conformally reseated metric g^v is physical. 



The original metric ^r^^ plays now the role of a variable providing an "already unified" fourth-order version 
of a theory including the gravitational metric plus a nonlinear scalar field (these theories form a basis 
for inflationary cosmological models). The vacuum Lagrangian (1.1) should first be transformed into the 
corresponding Einstein-frame Lagrangian (1.6), and then one has to add the interaction lagrangian for 
matter, minimally coupled to g^i,: 



L = 



i?-5^''0,^</.,,-2y(</.)+2Cat(«',ff) . (2.9) 



The matter Lagrangian ^mat is ^-independent since there is no physical motivation to assume that matter 
interacts with the scalar. Let us stress that in both cases ^mat is constructed in the same way: first one 
establishes by physical considerations the form of the free Lagrangian (or of the interaction Lagrangian for 
a number of coupled fields, e.g. for a complex scalar field minimally coupled to clcctromagnetism) for ^E" in 
Minkowski space and then one replaces the flat metric by the metric which is viewed as physical. Therefore 
•^matC^'j .9) in Case I and ^mat(^,.9) of Case II are the same fuctions of their respective arguments; clearly, 
for conformally-related metrics, £mat(^',.9) 7^ ^mat(^',5)- 

The field equations resulting from (2.9) now read 

= V(<^>5)+ V(*.5) (2.10) 

4-(^-atV^)=0 (2.12) 



6^ 

where f^ui/ is given in (2.5), with the potential as in (1.7) and (2.7). 

Here Tj^y = — -^7=1- 5^ (^mat a/— 5 ) > and due to the absence of any interaction between (f) and not only 



-g 

the total stress-energy tensor t^^+T^i, is conserved, but the stress tensors of each of the fields are separately 

conserved, V'^t^,,. = ^"T^u. = 0. 

To find the interaction of matter with the original metric g^i, one might naively trasform back the 

equation (2.10) to the Jordan frame, eliminating the scalar field cj) by the relations e^'^ = p = f'{R) and 
r\p{^)] = R{g); one would get 

f\R)R^. - i/(i?)<7M. - V;,V./'(i?) + □ f{R) = f\R)T^.{^, f{R)g), (2.13) 

where T^^ is the matter stress tensor defined in the Einstein frame and expressed in terms of the original 
unphysical metric g^i,. In general, however, equation (2.13) is incorrect. To see this, let us first recast the 
field equations (2.10) and (2.11) in terms of the variables {p,gixv), without assuming that r{p) = R{g). These 
read 

Gm.(5) =P"'(V^V.p-5^.np)-py(p)5M. + V(*,P5) , (2.14) 
„ 2p^dV 

Taking the trace of (2.14), eliminating □ p with the aid of (2.15) and using the explicit form of the potential, 
V = ^r{p) — 2^/[r(p)], one arrives at the equation 

R{9) - r{p)+g^''T^,{^,pg) = . (2.16) 

Recall that r{p) is defined by solving the equation f'[r{p)\ = p. In the absence of matter or if its stress tensor 
is traceless, (2.16) yields R{g) = r{p), then p = f'{R) and equations (2.13) do hold. In general, however. 
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setting p = f'{R) is inconsistent with (2.16). Provided that T^^ does not contain covariant derivatives of 5*, 
(2.16) can be viewed as an algebraic equation for the scalar field p. Solving this equation for p provides, in 
the presence of matter, the correct relation p = P{R; g, ^') which allows one to re-express the scalar field p 
in terms of the curvature scalar R{g) and obtain higher order equations of motion in terms of g and "I" only. 
To avoid confusion, let us call from now on the original unphysical frame, in which the vacuum Lagrangian 
takes the form L = f{R)^/--g, the vacuum Jordan conformal frame (VJCF), and the (also unphysical) frame 
into which the scalar field can be re absorbed, in the presence of the matter term in (2.9), the m,atter 
Jordan conformal frame (MJCF). Except for conformally invariant material systems, the two Jordan frames 
are different. A deeper understanding of the reason why the two frames do not coincide is provided by 
the Legendre transformation method. The explicit construction of the MJCF, the corresponding nonlinear 
Lagrangian and the resulting field equations are given in Appendix C. In MJCF, being the "already unified" 
frame, matter is nonminimally coupled to the metric and it is a generic feature that in the absence of 
gravitation, in flat spacetime, the nonlinear Lagrangian does not reduce to the standard form for a given 
species of matter. This fact significantly influences conservation laws. As is shown in Appendix C, it is 
possible to separate a "purely gravitational part" in the gravitational field equations in MJCF; this part 
satisfies the generalized Bianchi identity. Four matter conservation laws then follow from the equations in 
the same way as in General Relativity. These involve a number of terms mixing the curvature scalar with 
the matter variables and consequently do not resemble at all the elegant conservation laws T^i,'" = of 
Einstein's theory. One can only learn from them that in this frame the matter worldlines explicitly depend 
on curvature except for massless particles (photons), see Appendix B. 

In the Einstein frame, on the other hand, the Strong Equivalence Principle holds; the only possible 
difference between this version of gravity and General Relativity may arise from the physical interpretation 

given to the scalar field, whose entire role is confined to influence the metric field. In fact, (j) is assumed to 
describe a non-geometric spin-zero component of gravitation [22, 23]. 

2.3 Frames and conservation laws in STG theories 



Next we proceed to scalar-tensor theories of gravity (sec [1, 24] for recent reviews). These arc conceptu- 
ally different from NLG theories because these are not purely metric gravitational models, as the gravitational 
field is a doublet consisting of a spin-two field and a (non-geometric) spin-zero field, the Brans-Dicke scalar. 

The action of a generic STG theory is (we use conventions of Will's book [25]) 



S 



1 

16^ 



(fiR^ 



2(pA((^) 



t(*..9) 



(2.17) 



The "cosmological function" A(<^) is often omitted, and accordingly we consider here the models in which 
X{ip) = 0. If the coupling function u}{(f) is constant, the action is that of the Thiry-Jordan-Fierz-Brans- 
Dicke theory (see [1] and [25] for references). The field variables g^j/ and (p form the Jordan conformal 
frame (it is here that the concept originated). By assumption, ordinary matter minimally couples to the 
metric g^j, (and does not couple to the scalar gravity), thus by this assumption the Jordan frame is physical, 
i.e. describes measurable spacetime intervals. Proceeding as in the case of NLG theories one introduces a 
scalar variable 

1 1 

p=^=—L{g,^,^) = —ip; (2.18) 



-g OR ^' 167r' 

in this case it is not a function of the curvature but it coincides (up to a constant factor) with the already 
existing scalar field. With the aid of p one defines a new conformally-related metric ^j^i/ = PSfj-vt the 
transformation being already known since 1962 as Dicke transformation [26]. In terms of the new variables 
{Qhv-iV) the action is (up to a full divergence term) 



R 



S = j d^x^g 
and after a redefinition of the Brans-Dicke scalar. 



a;(</.) + - 



d(j) = uj{(j)) + 



(2.19) 



(2.20) 
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it takes the standard form of the action for the Unear massless scalar field minimally coupled to the metric 



S= (fx 



R- 



/ 167r Y 



167r _ 



(2.21) 



In terms of the Einstein-frame variables (g^i/,'/') the matter part of the action describes an interaction 
between ordinary matter and the scalar gravity. Clearly, it is due to the use of the unphysical (by assumption) 
variables for gravity. As a consequence the variational (with respect to ^^y) "energy-momentum" tensor 
for that can be defined by the Lagrangian ^^^imati'^, f^^gtiu) is different from the matter stress tensor 
defined in the Jordan frame and transformed to the Einstein frame; as already mentioned this notion is 
rather ambiguous and of little use. 

To study conservation laws for a STG theory one should first correctly identify the energy-momentum 
tensor for the spin- zero gravity. According to general rules of relativistic field theory it is provided by the 
variational derivative of the appropriate term in the Lagrangian. The latter should be identified with care 
since the interaction is nonminimal. First, the purely metric gravitational Lagrangian Lg should be separated 
out. To this end one formally views the Brans-Dicke scalar as a test field in a given fixed background 5^,^. 
The basic assumption of scalar tensor gravity theories is that the average value of the field (p determines the 
present value of the gravitational constant, < (p >= ^. Therefore, the situation where the scalar gravity 
is "switched off" does not correspond to = 0, but rather to assuming that <f is constant and equal to 
the present value of ^ (= Stt in our units); notice that iy9 = Stt is actually a solution of the field equation 
only if i? = 0. Setting = 8n and imat = reduces the Lagrangian in (2.17) to its purely metric part, 
Lg = ^R^/^ . This is the Einstein-Hilbert Lagrangian of General Relativity, as it should be expected. The 
scalar- gravity Lagrangian is thus defined (in the absence of matter) as, = L — Lg, i.e. 



L,„ — 



167r 



(^_87r)i?-^5'^V,^^,, 



(2.22) 



In the presence of minimally coupled matter the full Lagrangian in (2.17) is thus decomposed as L = 
Lg + L^+ (■mat^/'-g ■ The term L^ generates the variational energy-momentum tensor 



and the equation of motion 



. 



(2.23) 



(2.24) 



The invariancc of the action integrals Sm and S^p under spacctimc translations generates eight conservation 
laws (Noether's theorem, see Appendix A), V^t^^ = and VT^^ = (the latter can be directly verified 
using (2.24)). r^j^ explicitly depends on curvature. In equations (2.23) and (2.24) g^i, is an external metric 
field. When the back-reaction of p on the metric is accounted for, an ambiguity arises. Varying the action 
(2.17) with respect to g^i, one arrives at field equations for spin-two gravity. 



G^^{g) = Ti,^{ip,g)+T^^{^,g) , 



(2.25) 



and these allow one to eliminate the Einstein tensor from the expression for t^i, and the curvature scalar 
from (2.24). After these eliminations the field equations take the standard form which is usually applied in 
SGT theories. 



1 



Stt, 



Stt 

e^,{^p,g) + —T^,{^,g) , 
1 /- „ „ (hj ^ 



2a; + 3 



dip 



(2.26) 
(2.27) 
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By the basic assumption the scalar gravity docs not couple to ordinary matter and its role is c;onfined 
to influencing the metric, whereas (2.27) shows that the field ^p is influenced by matter although it does 
not cause a back-reaction. Actually interacts only with the metric field as is seen from (2.24) and 
(2.25). Furthermore, equation (2.26) defines the effective stress energy tensor 0^,^ for ip, acting as one of 
two sources for the metric field; this tensor is curvature-independent. Clearly, also t^^ is a source of metric 
gravity, according to (2.25). Thus, two stress-energy tensors are assigned to scalar gravity. Such ambiguity 
arises whenever spin- -two gravity is generated by two different sources since Einstein field equations alone 
determine only the total stress tensor, which can be expanded in various ways into contributions from each 
source separately. The effective stress tensor is not conserved and the relation 



yields 



V = V- + - 7) V (2-28) 

r-., = ^T^V,, . (2.29) 

It is worth stressing that 6^1, appears when the Brans-Dicke scalar acts as a source of metric gravity, while the 
variational definition of r^,,^ is always valid and the latter tensor should be viewed as the correct expression 
for a conserved stress- energy tensor (any possible ambiguities in the construction of are irrelevant in this 
respect). Whether or not r^^ provides a good physical notion of energy is a separate problem which will be 
addressed in Section 3. 

In studying conservation laws for matter in Einstein frame we restrict ourselves to the case of TJBD 
theory for simplicity, i.e. we set ijJ = const. Transforming (2.26) and (2.27) to the Einstein frame one finds 

the following field equations for g^^ and for the redefined scalar (f) = :^ ln(p, with 7 = (t^ + |) 

Gy^u = <t>,^<t>,u - \gf,ur'^(l>,c.<t>,i3 + 87re-^^T^,(*,5) 

= Vl*/-' 9) + 87re-^^T^,(*, g) , (2.30) 
□ <^ = A'K^e-'^'^T^^r" ■ (2.31) 

The notation T^,^(^',£() recalls that the matter stress- energy tensor is defined in the Jordan frame. By 
taking the divergence of (2.30) and using (2.31) one arrives at four matter conservation laws, 

V^T^. - 7r''(T^a0,/3 - = (2.32) 

Clearly the massless particles like photons or neutrinos follow null geodesic worldlines in both frames, while 
dust moves along timelike geodesic paths only in the Jordan frame (Appendix B). 

Finally, for the sake of completeness^, we comment on the conformally invariant scalar field model [27]. 
This theory fits to the general framework of dynamical systems generated from a general reativistic scalar 
field by means of conformal mappings, and is formally equivalent to a STG theory via redefinition of the 
scalar. However, it differs from STG theories in the physical interpretation given to the scalar field. The 
latter is in fact commonly viewed as a special kind of matter rather than being a spin zero component 
of the gravitational interaction. That the conformally invariant^" scalar field is equivalent to the massless 
linear field under a suitable conformal map is surprisingly little known to relativists although this fact was 
discovered by Bekenstein [28] twenty years ago^^. The form of the full action in Jordan frame. 



S = I cF^xyZ—g 



- ^5^"x,^x,. - ^x^R + Cat(*, g) 



(2.33) 



^ Actually one can also consider multiscalar tensor theories of gravity, see [1] 

^° The conformal invariance actually means the form-invariance of the field equation □ x~ ^Rx = under 
arbitrary conformal map 5^1/ 1-^ ^"^g^iu associated with x 1— > 0~^x (i-^- X is scaled like a particle mass). 
Only recently has appeared a work [29] applying the theorem. 
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clearly shows that this frame is assumed to be physical. It is natural to identify the Lagrangians for the 
metric and the scalar as Lg = ^^/^ R and = —^/—g{\g^''x,iJ.X,v + tlX^-^)) then the latter generates 
the variational stress-energy tensor for the scalar, t^i/(x,5). Instead of deriving it and the field equations 
from the Lagrangian, one views (2.33) as a version of STG theory. Upon comparing (2.33) with (2.17) and 
identifying = Stt — ^x^) finds 

Then the explicit form of T/^vix^Q) and the field equations follow from (2.23)-(2.25). As previously, the 

Einstein tensor and the curvature tensor can be eliminated, giving rise to the "effective stress-energy" tensor 
9fj,i, and to field equations analogous to (2.26) and (2.27). In the Jordan frame the Noether theorem implies 
T^^'^-i, = t'^"-v = 0, while (2.28) and (2.29) show that the effective stress-energy tensor 6^^ is not conserved 
in the presence of matter 

= -12x(6 - x^r^T'^^'x,. , X' < 6 (2.35) 

The condition x^ < 6 implies < </? < Stt and for these ranges of the field variables the theories (2.17) and 
(2.33) are equivalent to each other and to (2.21), since the conformal factor (2.18) is positive. The field x 

is redefined by d(/> = (l — gx^) ^ dx, '^-e. (}) — ^^In ^T^- '^'^^ inverse transformation back to the Jordan 
frame is 

= 2cosh^ ^-^^ and x = V^tanh ^-^^ (2.36) 

[For x^ > 6 ((/? < 0) the transformation from the Jordan to the Einstein frame is given by g^i, = —^^g^v 



and (f) — y|ln-^^-j-^; see Appendix E]. Thus the conformally invariant field x is nothing but yet another 

conformal image of the self-gravitating massless linear scalar (p. In general, a conformal map g^^ = F{<p)g^^ 
with arbitrary nonvanishing F transforms the Einstein-Hilbert Lagrangian for the scalar (j) minimally coupled 
to g^„ into 



(2.37) 



which can be transformed by further conformal maps and field redefinitions into any version of STG theories 
(Alonso et al. in ref. [11]). 

The interrelations between the various gravity theories are depicted in Figure 1. 

Wc (;onclude this section by emphasizing that the gravitational field equations for NLG in the Jordan 
conformal frame, (2.2) and (2.9), differ from each other in the matter part, and similarly equations (2.3) 
and (2.8) in the Einstein frame have different matter source term. Thus it is clear that matter dynamics 
depends on whether the physical metric is the one forming either the Jordan or the Einstein (or in any other) 
frame. Transforming from the assumed physical frame to the other (unphysical) one results in a number of 
bizarre terms depending on ordinary matter and/or the spin-0 gravity; nevertheless in each case one finds 
consistent conservation laws. In particular we stress that the divergence of the total energy momentum 
tensor (the sum of all the terms in the gravitational field equations depending on fields different from the 
metric one) does not provide a criterion for establishing which metric is physically acceptable. Accordingly, 
in the next section we will revert to a vacuum NLG theory (no ordinary matter), and study the distinguished 
role played by energy in gravitational physics; this will provide a motivation for regarding the Jordan metric 
as unphysical. 



Madsen [30] makes an incorrect statement on the subject. 
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3. Field redefinitions, physical variables and positivity of energy 

Having shown that, from the formal viewpoint, NLG and STG theories can be formulated in either of 

the two frames without giving rise to inconsistencies, and having clarified the theoretical implications of the 
choice of the physical metric, let us revert to our original question: can one choose the physical metric on 
arbitrary grounds? 

Brans' response [18] is negative: JBD theory is not "nothing but" General Relativity plus a scalar, and 

the metric tensor in the Jordan conformal frame has a direct operational meaning - test particles move on 
geodesic worldlines in this geometry. This is, however, a free assumtion and we have seen in Section 2 that 
there is no formal method which could determine which frame is physical. The problem is further obscured 
by the fact that there is no experimental evidence on interaction of the scalar with known matter. It might 
therefore seem (and implicit suggestions are sometimes heard) that a self-gravitating scalar field can be 
arbitrarily coupled to the spacetime metric. Presumably this is the origin of the view that physics cannot 
distinguish between conformal frames [16], the mere fact that the conformal mapping does not affect the 
particle mass ratios being clearly insufficient for proving it. 

Let us make the terminology more precise. Different formiilations of a theory in different variables 
(frames) will be referred to as various versions of the same theory. This includes not only mere transforma- 
tions of variables, but also transitions to dynamically equivalent frames. A theory (expressed in any version) 
is physical if there exists a maximally symmetric ground state solution which is classically stable. Classical 
stability means that the ground state solution is stable against small oscillations - there are no growing 
perturbation modes with imaginary frequencies. A viable physical theory can be semiclassically unstable: 
the ground state solution is separated by a finite barrier from a more stable (i.e. lower-energy) state and 
can decay into it by semiclassical barrier penetration [31]. The ground state solution may not exist, e.g. in 
Liouville field theory [32] , but in gravitational physics the existence of the ground state solution (Minkowski 
or de Sitter space) hardly needs justification. 

In most versions of a physical theory it is difficult to establish whether the ground state solution is stable 
or not and to extract its physical contents. Field variables (i.e. frames) are physical if they are operationally 
measurable and if field fluctuations around the ground state solution, expressed in terms of these variables, 
have positive energy density. Since the energy density cannot be defined for the metric field, the definition is 
directly applied to all other fields; the metric tensor influences the energy density of any matter and thereby 
the physical metric is indirectly determined. 

Energy density is sensitive to transformations of variables, particularly to conformal mappings. In terms 
of unphysical variables the energy density is indefinite and although the total energy is formally conserved, 
it loses most of its practical use. The ground state solution has (by definition) total energy equal to zero 
and when the theory is formulated in the physical variables (the "physical version") the solution represents 
the minimum of energy. Thus, stability is closely related to positivity of energy and instead of searching for 
growing perturbation modes one can study the total energy for nearby solutions [33] . 

These definitions apply to a relativistic classical field theory (and not to Newton gravity) and arc satisfied 
by all known unquantized matter. "For reasons of stability we expect all reasonable (though not quantum!) 
Seld theories to have positive energy density, and we expect all (classical and quantum) held theories to have 
positive mass-energy" [34]. The Weak Energy Condition is violated in some quantum states [35] while for 
all unquantized matter the Dominant Energy Condition holds [36]. Whenever the condition is violated one 
obtains physically meaningless results [37]. 

To show how this postulats work, consider a "ghost" complex scalar field in Minkowski space, minimally 
coupled to electromagnetism: 

L = --i-F'^^F^, + D^i,(D^i,Y , D^^ = + ieA^i' (3.1) 

iDTT 

Mathematically the theory is acceptable in the sense that the Cauchy problem for the field equations 

d^F^^ - 87rRe(ieV'*F"'V) = , D^D^^^Ij = (3.2) 
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is well posed. Physically, however, the theory is untenable since the full energy-momentum tensor 

T^u = ^ (j^^aF," - ^g^..F"^F„0^ - i [D^^{D^r) + D^D^r) - 9^.D^i^{D^xl,y] (3.3) 

is indefinite. The candidate ground state solution, V = F^i, = 0, has total energy E = and is unstable 

since any other solution with Ffj,^ = and Im(i/)) = has negative energy. The "ground state" decays 
via a self-excitation process where energy is pumped out from the scalar to the electromagnetic field and 
radiated away to infinity. Such a system, which can emit an infinite amount of radiation (while total energy 
is conserved) is clearly unphysical. 

Wc shall show that such effects do not occur for NLG thcoric^s if there is exact equivalence between 
the Jordan and Einstein frames. On the other hand, if the equivalence breaks down for some Lagrangians 
or solutions, the ground state solution (Minkowski space) is likely to be unstable. Before doing so, a few 
remarks on the problem of the physical (in)equivalence of the frames and the notion of energy are in order. 

A simple example showing that mathematical equivalence needs not imply a physical one is provided by 
classical Hamiltonian mechanics, where the canonical transformation Pj = and = —pi, being a mere 
renaming of positions and momenta, clearly shows that from the matematical viewpoint particle positions 
and momenta do not differ substantially. On the other hand, to construct operationally H{q,p) for a given 
mechanical system one should clearly discriminate positions and momenta, and kinetic and potential energy. 
In this sense, physical positions and momenta are those that are used to determine, on empirical and/or 
theoretical grounds, a Hamiltonian for the given system. Once H has been constructed in terms of physical 
variables, one has the freedom of making arbitrary canonical transformations. Whether or not the physical 
variables are uniquely determined is a separate problem and depends upon the system under question. 

To determine the physical variables and the Hamiltonian (or Lagrangian) for the system, it is in general 
insufficient to study the system alone: one should take into account its actual and possible interactions with 
its surroundings. The greater variety of interactions, the greater confidence that the dynamical variables 
describing the system and its Hamiltonian are correctly dc;fined. One; usually pays less attention to this aspect 
since in theoretical physics the physical variables are already given from empirical data and form a starting 
point for theoretical consideration: one is then interested in finding out the largest group of transformations 
for a given system, rather than in restricting the class of allowable frames. 

The very possibility that the system can interact with external agents means that the theory describing 
it is "open", in the sense that the surrounding in not included in the Lagrangian. On the contrary, in a 
"closed" theory the system constitutes the whole universe and no external agent can make an experiment 
on it; in this situation, any set of variables describing the system is equally physical and mathematical 
equivalence of frames means the physical one [19]. In a closed theory, in fact, energy is merely a first integral 
of motion without the distinguished features it has in an open theory. Each accepted physical theory is open 
in this sense (there arc some implicit trends in quantum gravity to view it as a closed theory [38]) and for 
instance in classical electrodynamics and quantum mechanics one has no doubts which variables have direct 
physical meaning and which are merely a convenient mathematical tool for solving a particular problem. 

In the case of gravity, in order to identify the physical variables and to formulate a physical version of 

the theory, one should experimentally study gravitational interactions of various forms of matter: motion of 
light and of charges and neutral test particles, behaviour of clocks and rigid rods etc. The point is that the 
presently available empirical data are too scarce to this aim. 

A theoretical criterion for pure gravity or for a system consisting of gravity and a scalar field is provided 
by energy, owing to its unique status in theories of gravity. In no other theories is energy effectively a charge. 
In any theory of gravity (including string-generated ones) the ADM energy should provide a good notion of 
energy [39] and the Positive Energy Theorem (see e.g. [40]) should hold. 

For a NLG theory one should compare energetics in the Jordan and Einstein frames. There is no 
generally accepted definition of gravitational energy for a higher-derivative theory, and in search for the 
physical metric one should not compare (as is usually done) the fourth-order version of the theory (1.2) 
with the second-order one. It turns out that the Einstein frame does not provide the unique second-order 
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version of the theory, and it was shown in [41] that any NLG theory (as well as theories with Lagrangians 
depending on Ricci and Weyl curvatures) can be recast in a version revealing a formal equivalence with 
General Relativity without changing the metric. (Yet the conformal rescaling of the metric is necessary 
to have a version of the theory with Einstein-Hilbert Lagrangian (1.6)). This is accomplished with the 
aid of a Helmholtz Lagrangian [41, 42] by applying a Legendre transformation. The Helmholtz Lagrangian 
contains the original unrescaled metric (the Jordan frame), and to get second-order equations of motion it 
is necessary to introduce a new independent field variable, a scalar field p. Thus, the field variables are now 
{9nu,p), and to distinguish the frame from the original Jordan frame (consisting of g^i, alone) it is referred 
to as Helmholtz Jordan conformal frame (HJCF). 

It should be stressed that the procedure is not an ad hoc trick and the scalar is (up to possible re- 
definitions) a canonical momentum conjugated to the metric g^y, and represents an additional degree of 
freedom existing already in a NLG theory, in comparison to (vacuum) General Relativity. The HJCF exists 
iff f"{R) 7^ (the i?-regularity condition), i.e. if the theory is a truly nonlinear one. 

The formalism [41] is outlined in Appendix C. The Helmholtz Lagrangian, which is dynamically equiv- 
alent to (1.1), is 

L„ = p[R{g) ~ r{p)]^ + f[r{p)]^ . (3.4) 

where, as previously, r{p) is a solution of the equation f'[r{p)] = p. The resulting field equations (C.3), after 
some manipulations, take on the form 

G^. = p-^V^V.p - i {p-'f[r{p)] + rip)} g^, = e^^{g,p) (3.5) 

and 

Dp=lf[r{p)]-^p-r{p) ; (3.6) 

these are equivalent to (2.14-15) in the absence of matter. Equations (3.5) have the form of Einstein field 
equations with the effective energy-momentum tensor of the scalar as a source^^. Although the scalar 
is a new independent variable, the number of degrees of freedom remains unchanged: the Lagrangian (1.1) 
describes a system with 3 degrees of freedom [22, 23] while it is obvious from the field equations (3.5) and 
(3.6) that Ljj represents two degrees of freedom for g^i, and one for p. 

Assuming that the Lagrangian f{R)y/^ does not contain the cosmological term, /(O) = 0, one finds 
that Minkowski space is a possible candidate groimd state solution to (1.2)^^. Then = ry^j^ and p = po 
with r{po) = -R(^) = are the solution in the HJCF. Is this solution stable? To assess it one considers 
solutions [g^iviP) to (3.5-6) which approach {rj^i,,Po) at spatial infinity at sufficient rate. Then the total 
energy of these solutions is given by the ADM surface integral at spatial infinity and the candidate ground 
state solution is stable (classically and semiclassically) if Eaum > and vanishes only for {rif^,^,po). In 
General Relativity the Positive Energy Theorem holds provided the r.h.s. of Einstein equations satisfies 
the Dominant Energy Condition for any source of gravity. The effective stress-energy tensor 6fj_i, for the 
spin-zero gravity does not satisfy the condition. The kinetic part of O^i, is quite bizarre: it contains second- 
order derivatives and is a homogeneous function of order zero in the field. The latter peculiarity can be 
removed by a field redefinition while the former cannot. In fact, setting p — F((p), with F' > 0, one has 
V^Vj,p = F'V^Vj,(/3 + F" dfj_ipdiyip and the term survives for any F. By setting p = e*^ the effective stress 
tensor becomes 

0^..{9, ^) = V^V.^ + d^ifd^if - ^ {e-'^/[r(e^)] + r{e^)} g^^ (3.7) 

The presence of a linear term in the effective energy momentum tensor for any long-range field is undesired 
b(X'aus(^ it causers difli( ulti(^s in determining the total energy. The general procedure to couple a field to 

The fact that an effective matter source arises from the nonlinear part of the Lagrangian was already 
noticed and applied in [48] in the case of a generic quadratic Lagrangian. 

We do not consider in this paper other possible ground states (for A 7^ 0), such as de Sitter and anti-de 
Sitter spaces, although our arguments can be suitably extended to deal with them. 
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gravity consists in regarding it at first as a tost fickl on a fixc;(l spacctimc background (which wc choose, 
for simphcity, to be Ricci-flat), whereby the energy-momentum tensor can be gauged by any identically 
conserved tensor a^v In the case of the scalar ip, the additional term ct^^ = V^Vj/t/? — g^i/D ip influences 
the total energy as measured at infinity. To avoid such ambiguity it is generally accepted that the energy- 
momentum tensor for any field should be quadratic in the field derivatives and contain no linear terms. In 
the case of equation (3.5) this argument is not crucial, because ^^j/ is uniquely determined by the original 
nonlinear Lagrangian and no ambiguity can arise. A definitive argument to regard the linear term in 
as unphysical is provided by the fact that due to its presence the energy density can be of either sign and 
be transported faster than light. The appearance of such linear terms signals that either the theory is 
unphysical at all (as is the case of Einstein- Gauss- Bonnet theory with compactified extra dimensions [43]) 
or merely that the field variables are unphysical and one should transfer to another frame. The same term 
also flaws the energy- momentum tensors (both variational and effective ones) for the Brans-Dicke scalar 
and the conformally-invariant scalar field in the Jordan frame. 

In terms of the HJCF field variables one cannot prove that solutions near the flat space have positive 
energy; on the contrary, the indeflniteness of O^i, deceptively suggests that negative-energy solutions exist 
and Minkowski space is classically unstable. To establish whether it is actually stable or not one makes the 
transformation to the Einstein frame. While the HJCF exists under the i?,-regularity condition f"{R) ^ 0, 
the Einstein frame exists if the further condition f'{R) > is satisfied at all spacetime points for a given 
solution (see Appendix E). Since one is interested in solutions which are asymptotically flat at infinity and 
do not differ too much (in the sense given below) from the flat solution in the interior, one requires f'{R) > 
and f"{R) ^ for i? ^ 0. These conditions hold for 

f{R) = R + aR'^ + f^kR^ . with a ^ . (3.8) 

fe=3 

The existence of the flrst two terms in the expansion, R + aR^, is of crucial importance for the equivalence 
of the Jordan and Einstein frames. These ensure the invertibility of the Lcgcndre transformation and yield 
p « 1 in the vicinity of flat space. The vicinity consists of all spacetimes for which R is close to zero; these 
include spaces of arbitrarily large Riemann curvature if e.g. ii^y = 0. The quadratic Lagrangian therefore 
carries most features of any NLG theory for which the equivalence holds. On the other hand, if (3.8) does 
not hold, the background solution - with respect to which the ADM energy is defined - has no counterpart 
in the Einstein frame, and therefore there is no hope that information on the total energy can be obtained 
by ordinary general-relativistic techniques. 

Let us assume that (3.8) holds. Consider now a spacelike 3- surface S embedded in an asympotically 
flat spacetime (Af, f;,;,^). The asymptotic flatness in the HJCF is defined as in General Relativity since the 
field equations are of second order; the weakest assumptions [44] are (in obvious notation) 

9^,. = ri^u + 0(r-i-^) , = 0(r-t-^) with e > . (3.9) 

Then the leading-order contribution to is i? = 0{r~i~^). The HJCF dynamics implies p = f'{R) and the 
Lagrangian (3.8) yields the asymptotic behaviour of the scalar field, 

J9= l + 0(i?) = l + 0(r-5-<^) . (3.10) 

Under these conditions one now proves, contrary to what might be expected from the properties of 9^^, that 
the total energy is positive for a NLG theory. 

Positive Energy Theorem for NLG Theories. Let S be an asymptotically Bat, nonsigular spacelike 
hypersurface in a spacetime (M, gr^,^) topologically equivalent to JR*. If: 

i) the Lagrangian is given by (3.8); 

a) p> and f"{R) ^ everywhere on E; 

Hi) a solution {gfj,v,p) to (3.5-6) satisfies the condition (3.9) (and hence (3.10)) on S; 
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iv) the potential V (1.7) for the scalar </> in the Einstein frame is nonnegative, V{(t)) > 0, on S; 
then the ADM energy in the Jordan frame is nonnegative, 

EAuuig] = ^ / dS^ig.j^j - gjj^i) > . (3.11) 

Proof. The proof is a direct extension of that given by Strominger [20] for / = i? + aE?. Consider the 

total energy of the conformally related solution {gp,u,4') iii the Einstein frame. Because of the faUoff rate 
in the HJCF, the metric g^^ = pg^u is asymptoticahy flat and the total energy of the solution is finite and 
given by the ADM integral over a boundary 2-surface at infinity, EAuuig] = \ jgi dSi{gijj — gjj^i). If 
V > the Dominant Energy Condition holds and this energy is nonnegative [40]. Now, replacing by gf^^, 
and applying (3.10) one easily finds that EABuig] = -£'Adm[5] and hence that i^ADiviis] > 0. □ 

There is a subtle conceptual difference between our proof and Strominger's one: the use of the ADM 
integral for the total energy is not based on the fact that at large distances the dynamics is governed by the 
lower-derivative terms in (1.2), but rather on the existence of the HJCF. In the latter the integral is defined 

as in General Relativity. 

It should be stressed that the global equivalence between the Jordan and Einstein conformal frames 
holds if p > and f"{R) ^ everywhere on (M, g^^), while to prove the theorem one needs only to assume 
that these conditions hold on the initial data surface S, what implies that the frames are equivalent in some 
neighbourhood of S. Consider for example the collapse of a cloud of dust. The exact relation between R 
and the dust energy density p depends on the form of the field equations (which metric is physical. Section 
2); in general one expects that R grows when p does. On the initial surface E the dust is diluted and 
i? ~ 0, thereby the frames are equivalent in the vicinity of S. Near the singularity p is divergent and R 
is unbounded, p may change sign and thus in this region of the manifold there may be no mapping to the 
Einstein frame. Yet the energy, being conserved, is still given by its value on S. 

In the class of solutions for which the two conformal frames are at least locally equivalent, Minkowski 
space is the unique one having zero energy. In fact, let M contain a spacelike surface S on which the 
assumptions of the theorem hold. Then, they hold in some neighbourhood ZV(S) and the spacetime region 
{U{T,), g^i,) can be mapped onto (W(E), gp,y). Let -BadmI,?] = when evaluated on S; then g^,^ is flat 
in U{Ti) (Minkowski spacetime is the unique zero-energy solution in General Relativity) and the solution 
{Q/j-v = Viii^j <!) = 0) is identically mapped back onto (g^i, = = 1) in U{T,). Then the spacetime (M, g^^) 

is flat in the open region W(E) and this solution can be analytically extended onto the entire manifold M, 
thereby the spacetime actually is Minkowski space topology is always assumed). 

Thus the total energy is positive in the Jordan frame for all solutions to (1.2) containing an open region 
U{±) where these are close to Minkowski space {R close to zero) and where y(^) > 0. The flat solution is 
classically stable against decay into these solutions and is a ground state solution for a given NLG theory. 

Now wc consider a few examples. 

(1) f{R) = R + aR^. The Lagrangian is R- regular, since /" = 2a 7^ everywhere. 

(A) a > 0. This is the case studied by Strominger [20]. Using the variable p instead of (j) one has 
V{p) = 8^(P — 1)^ > everywhere, while p > for i? > — For all solutions in the Jordan 
frame such that there exists a spacelike asymptotically flat surface E with R> — ^ on it, the total 
energy is nonnegative and flat space is classically stable. Strominger has also shown that Minkowski 
space is the unique solution in the Jordan frame with vanishing energy. 

(B) a < 0. The conformal equivalence holds in the region R < \^\, but V{p) < everywhere and the 
theorem does not hold. One may expect that in the vicinity of flat space there are solutions with 
negative energy and Minkowski space is classically unstable. 

(2) /(i?) = ^(6"^ — 1). This is again a i?-regular Lagrangian. Furthermore, the Jordan and Einstein frames 
are globally equivalent since p = e"^ > for any solution. The potential V{p) = 2^(lnp + ^ — 1) is 
positive for a > and negative for a < {V{1) = 0). Thus for a > all asymptotically flat solutions 
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have positive total energy and flat space is stable against any perturbations, small or large (classical 
and semiclassical stability) . On the contrary, for a < the Dominant Energy Condition does not hold, 
what signals existence of negative-energy states and classical instability of flat space. 

(3) f{R) = i ln(l + aR). In this case one has f"{R) = -jy^^ 7^ for i? ^ -i, and V{p) = ^{Inp - 
p+1). For a > one finds V{p) < for p ^ 1 and flat space can classically decay into negative-energy 
states. For a < the conformal equivalence holds for i? < — i and V{p) > for all p ^ 1, thus for these 
solutions the total energy is positive and flat space is classically stable. 

In all these examples states close to Minkowski space have positive energy if the coefficient of the R^ term 
in the Taylor expansion of / is positive. This is a generic feature of Lagrangians (3.8). For R close to 
zero the two frames are equivalent and the theorem holds iff the potential is positive. For vacuum the 
relation r{p) = R{g) holds; using it in the expression for V{p), given after equation (2.15), and expanding 
the potential aroimd R = one arrives a.t V = ^aR^ + 0{R^) for solutions. All R = solutions represent a 
local extremum of the potential. Hence energetics and stability of ffat space are determined by the i?^ term 
in the Lagrangian of a NLG theory. 

The lowest-order contribution R+aR? to the full Lagrangian is crucial for classical stability of Minkowski 
space. Whenever these terms are absent, e.g. for / = R^ , fc > 2, the conformal equivalence with the Einstein 
frame version is broken at flat space. In this case as well as when the Einstein frame is defined for i? = 
but the potential is negative on S, the classical decay of flat space may occur because the Positive Energy 
Theorem at Null Infinity (see [40]) does not hold. While the ADM energy is conserved, there may be metric 
perturbations about the background such that gravitational waves carry away unbound amounts of energy 
and the Bondi-Sachs mass remaining in the system decreases to minus infinity. 

Whenever the ADM mass is positive for solutions close to flat space (i.e. if a > 0), the latter need not 
be semiclassically stable. There may exist solutions (in the Jordan frame) with negative energy which are 
far from Minkowski space. These are spacetimes so highly curved (large \R\) in the interior, that either the 
equivalence with the Einstein frame is violated {p changes sign on each asymptotically flat S) or the potential 
becomes negative in a region of any S. The occurrence of any of the two possibilities depends on the form 
of the Lagrangian, e.g. for the exponential one (example 2 above) there are no such solutions, while for a 
cubic Lagrangian, R + aR^ — the potential is negative for large p and some values of a. Semiclassical 

instability of standard ground state solutions is a generic feature of higher-dimensional theories: it was 
found both in Kaluza-Klein theory in five [31, 45] and in ten [46] dimensions, as well as in string theory 
[46]. Although disturbing in itself, it is not dangerous. Classical instability is a more severe problem and 
usually makes the theory untenable. If there is another possible ground state which is stable, however, it is 
harmless, as in the case of the uncompactified Einstein-Gauss-Bonnet theory, where rf-dimensional anti-de 
Sitter space is unstable while Minkowski space is stable [47]. Minkowski space is classically unstable in a 
gravity theory with Lagrangian aR + bR^ + cCa/Sfj.i^C"'^'^'^ [48] and in a four dimensional theory resulting 
from Einstein- -Gauss-Bonnet theory with compactified extra dimensions [43] (in all these cases there is no 
conformal mapping of the theory onto a model in General Relativity). 

Equality of the total energy in the two conformal frames is due to the fact that -Eadm , being effectively 
a charge, is evaluated at spatial infinity (where p 1) and is rather loosely related to the interior. The only 
detailed information about the interior of a gravitating system that is needed is whether all local matter 
energy-flow vectors are timelike or null. This connection is lost in the Jordan frame. Also the Bondi-Sachs 
mass, when expressed in terms of the Jordan frame, does not satisfy the correct evolution equations, as is 
mentioned in [1]. These flaws do not reflect the genuine properties of spin-0 and spin-2 gravity, these are 
merely due to an improper choice of field variables. The effective stress tensor 6*^1/ (3.7) does not provide 
the true physical energy-momentum of the scalar field. As any redefinition of the scalar cannot help, it is 
the metric variable which needs redefinition. 

Accordingly, the Jordan frame should be regarded as unphysical not by an arbitrary choice of definition 
but because in the second order version of the theory - whereby the spin-0 and spin 2 component of gravity 
are singled out - these variables are unrelated to the total energy of the system. On the other hand, the 
Einstein frame satisfies all general requirements of relativistic field theories. Therefore it is the physical 
frame, and g^,^ determines the spacetime intervals in the real world. It provides the physical version of a 
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NLG theory. All other conformally related versions of the theory, including the initial one, are unphysical in 
some aspects^^. This physical frame is uniquely determined by the Lagrangian (1.1) and the physical metric 
for the vacuum theory is = f'{R)g^v Physically, this transformation means a transition to the frame 
where all fields, including spin-0 gravity, satisfy the Dominant Energy Condition. Then it is possible to 
establish that Minkowski space is a classically stable local minimum of the energy. This holds true, however, 
if the potential of the scalar field does not attain a local maximum at flat space, which would destabilize the 
latter. The only remedy is to exclude as unphysical all Lagrangians giving rise to negative potentials near 
flat space, and this amounts to the requirement a > in (3.8). 

4. Conclusions 

Our analysis of the physical features of energy in NLG theories leads to the conclusion that the Einstein- 
frame metric should be regarded as the physical one. In this perspective, whenever a nonlinear Lagrangian L 
admits flat space as a stable ground state solution, it is physically equivalent to a scalar field with a potential 
V (determined by L) minimally coupled to the rescaled metric g^i, in General Relativity. Then, any form 
of matter should minimally couple to g^i, and not to the scalar. The Strong Equivalence Principle is then 
retained and the scalar, which can (though not necessarily) be viewed as a spin-0 component of universal 
gravity, appears only as a contribution to the full source in the Einstein field equations. 

Scalar fields are now very popular in cosmology. Instead of introducing ad hoc a special form of the scalar 
field potential in order to solve a particular problem, one can start from some nonlinear Lagrangian which 
provides the desired potential, provided that one can give a deeper motivation for the former. Our conclusion 
does not mean that any viable nonlinear gravity theory is identical with General Relativity. Although the 
scalar field has no impact on motion of ordinary matter (in particular, free test particles move on geodesic 
lines), the space of solutions is different. For instance, for nonconstant scalar field there are no black holes 
with regular event horizons since a black hole cannot have a scalar hair [49]. 

It is amusing to notice that twenty years ago Bicknell [2] concluded his study of the phenomenology of 
purely quadratic theories with the words: "These results eliminate gravitational theories based on quadratic 
Lagrangians from the realm of viable gravitational theories and point strongly towards the uniqueness of 
the Einstein equations" . Our results regarding the Lagrangians giving rise to potentials which attain a 
maximum for flat space (a < 0) suggest the some conclusion. However, a rigorous proof of the instability of 
flat space in that case is not yet available. The case of Lagrangians for which a conformal mapping onto the 
Einstein frame does not exists at Minkowski space (e.g. i?'^) is more subtle. The presently known techniques 
of investigating energy and stability fail there and it remains an open problem whether positivity of energy 
can be shown for a subclass of these Lagrangians. Another open problem is how to deal with Lagrangians 
which are not of class C'^ at i? = 0, e.g. R + aR^ + i?3 . 

Finally, it should be emphasized that our method of determining the physical metric for a gravity theory, 
by studying the energy and stability in the vacuum theory, does not apply when the interaction of matter 
with gravity is already prescribed by a more fundamental theory. This is the case, e.g. of string-generated 
gravity, where the string action in the four-dimensional field-theory limit contains a number of fundamental 
bosonic fields coupled to the metric fiels and to the dilaton. The problem of identifying the physical metric 
appears there too [50] and should be dealt with in a separate way. 
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Appendix A: Generalized Bianchi Identity and Conservation Laws 



Although it is well known to the expert in gravitational physics, we provide here, for the reader's 

convenience, the derivation of the Nocther conservation laws for gravitating matter, based on a generalized 
Bianchi identity. Consider a generic action for gravity and a matter field 



S — Sg + Sm — I d^X\/^ F{gcef3, Rafjpa, ^ ixRafipa, ■ ■ ■ , ^ fj-i ■ ■ ■ ^iJ.„Raf3pa) 



(^.1) 



here -F is a scalar function depending on the metric and its derivatives (via the curvature tensor) up to 



(n + 2)th order. One takes an infinitesimal point transformation, x'' 



+ e^'^{x), with S,'^ being an 



arbitrary vector field which vanishes, together with all its derivatives up to (n + 2)th order, at the boundary 
of the region il. The action integrals for gravity and matter, Sg and Sm, are separately conserved under the 
transformation. The variation of the gravitational action is 

5Sg = = [ d^x5{^/^ F) . (A2) 

By applying (n + 2) times the Gauss theorem and discarding the surface integrals one arrives at the standard 
formula 



where 



SS„ 



I 

Jn 



da 



+ ••• + (- ir+'5„,---a„„+ 



Sg' 



{A.3) 



+ dadfi 



dg'"' 



,ai---a„-|-2 



One replaces the tensor densities by pure tensors according to 



then 



6Sa 



Jn 



{AA) 

(^.5) 
(^.6) 



For the case of the infinitesimal point transformation one finds that Sg'^'^ = e{^^''^ +^'^'^), what ensures that 
the surface integrals in 5Sg vanish. Furthermore, 

= 5Sg = e [ d^x^Q^^e'" + r") = 2e / d'x^Q^^^'" 
Jn Jn 

=2e / d^x^iW'iQ^^i^') - ^'^V^Q^.] = 2e / Q^.^d^^ - 2e / d^^cv^^Q^.'^ 
Jn Jdn Jn 

The surface integral vanishes, and due to the arbitrariness of the vector inside the integration domain fl 

one gets Qfw '^ = 0. This is an identity which holds for any metric, independently of the field equations. 
In the case F = f{R) the tensor Q^^ is equal to the l.h.s. of equation (1.2). Any scalar function F of the 
spacetime metric and its partial derivatives up to (n + 2)th order gives rise to a generalized Bianchi identity 



1 5{^F) 



. 



(^.7) 



The symmetric tensor Q^i, is a function of the metric and its derivatives up to (2n + 4)th order, unless 
F is a linear combination (with constant coefiacients) of particular expressions like R (Einstein-Hilbert 
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Lagrangian), RaPpaR"''^'''^ — '^RadR"'^ + R^ (Gauss Bonnet term), etc. (which correspond to the Euler 
Poincare topological densities in dimension two, four, and so on). In other terms, we can construct out of 
the Riemann tensor, besides the well-known contracted Bianchi identity G^y'" = 0, an infinite number of 
differential identities (which can all be derived from the ordinary Bianchi identity and its higher derivatives), 

sometimes called strong conservation laws. 

Upong variating the action (A.l) with respect to arbitrary variations 6g^'^ of the metric one arrives at 
the gravitational field equations 

The matter stress-energy tensor, defined by y/^T^^ = — 2^J^(y^^niat) may depend on first and second 
derivatives of the metric. The coefficient ^ in (A. 8) is due to the normalization of the gravitational Lagrangian 
assumed in (^.1). Actually, in order to have correspondence to General Relativity, upon expanding F in 
powers of the curvature scalar and other curvature invariants, if the linear term R has a coefficient a, then 
the standard matter Lagrangian ^mat in (^-1) should have coefficient 2a (having set SttG = 1). Assuming 
that the gravitational field equations hold, one takes the divergence of (A. 8), and using the identity {A.I) 
one gets four conservation laws for matter, T^^'"^ = 0. 

An equivalent and somewhat shorter way of derivation of the conservation laws consists in utilizing the 
invariance of the matter action Sm under infinitesimal point transformations. These transformations give 
rise to the variations Sg'^^ — 2eS}-'^'^^ and 5"^ . Assuming that the variation 8'^ and its relevant derivatives 
vanish at the boundary 9f2 of the integration domain, one finds that 

SSm = 0= [ d'^x 
Jn 

Provided that the equations of motion for matter hold, = 0, the invariance of the action yields 

SSm = -\ I d^xV^T^Jg^'' = -e f d^x^T^^^''" = e / d^xV^eT^.u" , (A.IO) 
^ Jn Jn Jn 

hence T^jy'" = 0. These conservation laws are actually either equivalent to the equations of motion for matter 
or to a subset of them (the number of equations of motion may exceed four). 



6^ 



{A.9) 



Appendix B: Motion of dust and photons in the Jordan and Einstein conformal frames 

We consider here the explicit forms of the conservation laws in the cases I and II for NLG theories 
and in the Einstein frame for STG theories. We show on explicit examples that the different versions of 
conservation laws reflect the known geometrical property that null geodesic worldlines are preserved under 
conformal rescaling, while timelike geodesic paths in one frame correspond to non-geodesic paths in the other 
frame: this fact has a physical counterpart in the different coupling between matter and scalar field in the 
two frames. The consistency of the conservation laws described in Section 2 with the physical assumptions 
is thus ensured both in case I and in case II. 



Case I: the Jordan frame is physical. 

In Einstein frame the conservation laws for matter are given by (2.8). Inserting the explicit form of the 
(jf)-field stress tensor t^v and making use of the equation of motion (2.6) one arrives at four equations 



gap 



^aT^p + —=Tap4>,ij, 

v6 



= 



(B.l) 



Clarly, as is seen from the Lagrangians (2.1) and (2.3), only conformally invariant matter does not interact 
with the scalar (for this matter g^^T^p = and the equation of motion for ^ becomes ^'-independent). This 
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is the case of electromagnetic field, hence photons travel along null geodesic lines both in the physical and 
in the rescaled geometries. We show that this result also follows from the conservation laws (B.l). 

Consider the null eletromagnetic field which is the curved- space classical representation of the photon, 

F^i, = k^q„ - k„q^ , with k^^k^ = k^^q^ = ; {B.2) 

here, A:^ are the covariant components of the wave vector and q/^ is a spacelike covector determined up to 
Qn ^ Qfi + ^kf^. Since — d^Ai, — di^A^ is metric-independent, the components k^ and q^ are independent 
of the conformal rescaling of the metric. Introducing in Jordan conformal frame a scalar w = -^g"^qaQi3 > 
one finds that the stress tensor for the null field is given in this frame by Tf^„{F, g) = w k^k^. Then WvT'^" = 
imply that k'^V^k^ oc k^ - the wave vector is tangent to null geodesic curves. 

In transforming to the Einstein frame, pg^u, k^ 1-^ k^ and i-> g^, thus Tn„{F,p ^g) = pwkfj,k„, 

with w = -^g"^qaQf3- Equations (B.l) for this stress-energy tensor reduce to 



k^y^k" 



and the wave vector remains tangent to null geodesies but the parameter a defined by k^ = is not the 
afJine parameter along them. 

On the other hand, massive particles do interact with the (/)-field. Consider a self-gravitating dust. In 

Jordan frame its stress tensor is Tp,^ = pu^u,^ and \7j^T^'^ = gives rise to motion along timelike geodesies, 
vy\l = 0. Under the conformal rescahng one finds that = p^u^ and p = p^p [28], thus the dust stress 
tensor defined in the Jordan frame and expressed in terms of the Einstein frame While in 

the Jordan frame the density current is conserved, 'Vv{pu'^) = 0, in Einstein frame the corresponding density 
current is coupled to the scalar. 

Using this fact in the conservation laws {B.l) one arrives at the following equations of motion for dust in 
g^i, metric: 

v6 

In particular, the dust worldlines are geodesic for only if = 0. 

Case II: the Einstein frame is physical. 

Consider a pressureless dust either being a source in the field equations (2.10) or just a cloud of test 
particles. In the Einstein frame its energy-momentum tensor is divergenceless and the particles follow 
timelike geodesic curves. In MJCF the interaction with gravity becomes highly complicated (equation ((7.19) 

in Appendix C), and the worldlines are not geodesic. 

Massless non-interacting particles move along null geodesic paths in the physical metric g^^, and this 
property is preserved under any conformal rescaling. This fact can be also explicitly derived from the matter 

conservation laws valid in MJCF. For conformally invariant matter-'^^ the field equations (C19) reduce to 
(2.13) and the generalized Bianchi identity gives rise to the following four conservation laws: 

V^[/'(i?)T^.(*, f{R)gc.f3)] = . (B.5) 

Thus the proof essentially follows the same lines as in Case I. 



In this case the MJCF coincides with the VJCF. 
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Scalai — tensor gravity theories 

By assumption, in the Jordan frame dust particles follow timelike geodesic worldlines and photons follow 
the null ones. In the Einstein frame, we can relate the motion of matter to the conservation laws (2.33) by 
the same procedure used for NLG theories. Case I; in fact, equation {B.l) coincides with (2.33) with w = 0. 
The calculations for the general case follow essentially the same lines, leading to the conclusion that photon 
worldlines are geodesic in both frames while dust particles move along nongeodesic curves for the metric 5^^: 

u^V^u^" = i7(<^''^ + 0,.w''w") . (B.6) 



Appendix C: Direct and Inverse Legendre Transformations in the Presence of Matter 

The equivalence of any NLG theory (1.1) to General Relativity plus the scalar field (1.6) is dynamical 
in the sense that the spaces of classical solutions for each theory are locally isomorphic (see Appendix E), 
the isomorphism being given by (1.3). This is directly achieved by the conformal transformation of the 
field equations (1.2) [2, 5], and is suSicient for various purposes. On the other hand, using this ad hoc 
procedure it is rather difficult to establish the particle contents of NLG theories [22, 23] and to recognize 
the field (j) as an independent degree of freedom present in the theory; instead one merely views as a 
function of the curvature scalar R of the original metric [6]. These goals can be easily achieved using a 
Lagrangian formalism. A Legendre transformation allows to transform the NLG Lagrangian (1.1) into a 
dynamically equivalent one, linear in the curvature scalar and including an auxiliary scalar field. This 
equivalent Lagrangian is named "Helmholtz Lagrangian" by analogy with classical mechanics. Besides 
being elegant and mathematically well-grounded, the Legendre transformation and the related Helmholtz 
Lagrangian turn out to be indispensable tools in dealing with NLG theories: one would hardly guess how 
matter fields with a prescribed coupling with the Einstein-frame metric interact with the Jordan-frame 
metric, without relying on the corresponding Lagrangians. As is shown in Section 2, a naive approach leads 
to inconsistencies. 

Being linear in the second derivatives of the metric tensor, the Helmholtz Lagrangian generates second- 
order gravitational equations; then the Einstein frame Lagrangian (1.6) is obtained by re- expressing the 
former in terms of the conformally rescaled metric g. It is natural then to use the Helmholtz Lagrangian to 
compare the dynamics in the two frames, and therefore it is a central ingredient in OTir approach. In this 
Appendix we outline of the general setting of the method [3, 4, 41] (we refer the reader to [51] for a rigorous 
mathematical exposition), and then present in detail how matter interaction terms should be inserted into 
the gravitational Lagrangians (both the original nonlinear one and the Helmholtz Lagrangian) when the 
matter fields are assumed to be minimally coupled to the Einstein-frame metric. 



C.l Helmholtz Lagrangian. 

For the nonlinear lagrangian density (1.1) one introduces the generalized conjugate momentum 



The scalar field p is clearly nothing else but the scalar defined by (1.3). As in Section 1, let r{p) be a function 
such that f ji^^.f^p-j = Pi such a function exists if f"{R) 7^ 0, and in this case we say that the nonlinear 
Lagrangian (1.1) is R-regular. The Helmholtz Lagrangian corresponding to (1.1) is then defined as 

= p[R{g) - r{p)]^ + Mp)]^. (C.2) 



The action functional corresponding to (C.2) is formally a degenerate case of the STG action (2.7), with 
w = and nontrivial cosmological function. A basic feature of the Helmholtz Lagrangian is that it does not 
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contain derivatives of the scalar field p. The Euler Lagrange equations obtained by varying the metric g^i, 
and the scalar p independently in the action defined by (C.2) are 

R{g) = r{p) 

(C.3) 

,pGtiv{g) = "^tNvP- gtiunp+\{f[r{p)]-p-r{p))g^^^ 

which are manifestly equivalent to (1.2). 

What has been done is nothing but a generalization of a standard method of classical mechanics. Given 

a first-order Lagrangian L = L{t,q'^,q^) for a system of point particles, the second-order Euler-Lagrange 
equations — §^ = can be recast into a first-order system, provided the Lagrangian is regular 

(det Q^iQqj 7^ 0). This can be done in two ways: either one writes the Euler-Lagrange equations in normal 

form, = a^{t,q^ ,q^), then one introduces independent velocity variables and writes the equivalent 
system in the velocity space 

ii^ = a'{t,q^,u^) ^'^'^^ 

or one defines the Legendre map {q^,u^) i-^ {<t,Pi = ^) and finds the inverse map = u^{q^,Pj), which 
transforms the previous first-order system into the Hamilton equations in the phase space of the system. An 
interesting point, which is not always emphasized in textbooks, is that both first-order systems arise from a 
variational principle, defined by a degenerate first-order Lagrangian, the Helmholtz Lagrangian [42] 

= p, {{i - u') + L{q\u^) , 

which can be regarded either as a function over the velocity space, whereby {q\u^) are the dynamical 
variables and pi = pi{q^,u^), or as a function over the phase space, where the variables are {q^ ,Pj) and 
— u^{q^ ,pj) is the inverse of the Legendre map defined above. In fact, the solutions of (C.4) are the 
extremal curves in the velocity space for the action functional 

S,= J {Piiq',u^) [q' - u'] + L{q\ u^)} dt , 

while the Hamilton equations can be derived from the action (in phase space) 

S^= J {pi [e-u\q^,pj)] +L[q\u^{q\pk)]}dt. 

The Lagrangian (C.2) plays the role of the Helmholtz Lagrangian in phase space for any i?-regular NLG 
Lagrangian (1.1). As in particle mechanics, it is also possible to define a Helmholtz Lagrangian in the velocity 
space: this corresponds to regarding directly the "velocity"-/? as an independent field u (as is done e.g. in [6], 
following [52] and without invoking the notion of Helmholtz Lagrangian), and introducing the Lagrangian: 

= f'{u)[Rig) - u]^ + /[«]^/^. (0.2') 

The corresponding field equations are nothing but the system (C.3), upon the substitution p f'{u), 
r{p) u. For quadratic Lagrangians, which most frequently occur in the literature, the two fields u and 
p coincide up to a constant factor. More generally, since the i?-regularity is still necessary to ensure the 
equivalence of (C.2') and (1.1), the choice of either p or u to represent the independent scalar field occurring 
in the second order picture is, from the mathematical viewpoint, a mere matter of convenience. The use 
of the scalar field u allows one to bypass the problem of finding explicitly the function r{p), but then the 
gravitational equation contains a rather complicated dynamical term for the u-field, V^V^/'(tt)— g^j/ □ f'{u), 
which depends on the particular choice of f{R) in (1.1). On the other hand, there is a physical motivation 
to formulate the theory in terms the scalar field <p, which can be equivalently defined either by cx Inp 
or by </) oc In f'{u), since this field interacts with the Einstein-frame metric in the standard way (to this 
purpose, the problem of inverting /' cannot be avoided). In the sequel of this Appendix, however, it will be 
convenient to write all the equations in terms of the conjugate momentum p. 
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C.2 Inverse Legendre Transformation. 

So far, we have recalled how a nonlinear Lagrangian can be recast into a Helmholtz Lagrangian. Now 
wc present the inverse procedure: how to find a nonlinear Lagrangian equivalent to a given Helmholtz 
Lagrangian. In fact, let 

L,{p;g,^) = ,/^\pR-H{p;g,^)] (C.5) 

be a generic Helmholtz Lagrangian. The function H plays the role of a Hamiltonian^'^ for a system of matter 
(denoted collectively by 5*) and gravity g^y, thus it does not depend on derivatives of the gravitational 
momentum p, while it depends on covariant derivatives of \1/ up to some order (the semicolon separates the 
field variables which are accompanied by their derivatives from those which are not). The variation of the 
action = J L^d^x with respect to ^ yields the equations of motion for matter, 

^{^gH)=Q, (C.6) 
while variation of it with respect to p yields an algebraic equation for the canonical momentum, 

Hi,)-'-S^ = o. (a7) 

Any solution of this "equation of motion" for p is denoted by P{R;g,^). Finally, varying the metric one 
gets (second-order) gravitational field equations 



1 (5 1 

— {yf^pR) =pG^,^{g)-Vij.VuP + gf,iyDp-- 



/=5 Sg^' 



-9 Sgi^" 



iV^H) ; 



(C.8) 



clearly (C.7) and (C8) are a generalization of (C.3). One now defines, with the aid of a solution to (C.7), a 
nonlinear Lagrangian for the metric and the matter: 



P{R;g,'^)RV^ -H[P{R;g,'^y,g,'^]V^ , 



(C.9) 



where R stands for R{g) (and is not and independent variable). To show the equivalence of the two La- 
grangians, one finds the field equations resulting from the stationarity of the action 



5nl = / L,Ld^x = J d^xV^[P{R;g,'i!)R-W{g,n , 
where W{g, ^) = H[P{R; g, g, For a variation of the matter fields, one has 



(C.IO) 
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^1 



m + 



n=l 



dp 



dP 



P=P J 
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(5* 



hi 



d^x 



(5* 



p=P 



hence the equations of motion for matter are 



This "Hamiltonian" has actually nothing to do with the notion of energy and the ADM formalism in 
General Relativity. 
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(5* 







(C.ll) 



p=p 



and are equivalent to (C.6), since (C.7) holds identically. The latter plays also a crucial role in deriving field 
equations for the metric. In fact, a variation Sg^" in the action yields 



= J d^x[Pd{y 



-gR) 



-gRSP-SiV^W)] 



and 



-g [R - — 



dp 



sf^R5P-5{^W) = 
so that the equations for gravity are 

P{R; g, ^)GM - "^y^^^P + g^uUP 



SP- 



5H 



Sg"" 



p=p 



1 



SH 



-9 Sgt'" 



{C.12) 



p=p 



and these are equivalent to (C.8). This shows that the Lagrangians and equivalent, at least 

whenever equation (C.7) has a unique solution (see Appendix E). 

One now uses the general formalism described above to find out a nonlinear Lagrangian, and the 

corresponding fourth-order field equations in MJCF in Case II of Section 2, i.e. when matter is minimally 
coupled to the rescaled metric g^i^. The inverse transformation from the Einstein frame does not lead back 
to the original VJCF but to another (conformally related) spacetime metric. This makes no surprise since, 
after adding interaction with matter in the Einstein frame, the inverse transformation is applied to a different 
system. Let us sec it in more detail. 

C. 3 Adding matter interaction in the Einstein frame. 

The starting point is a nonlinear vacuum Lagrangian for pure gravity. Denoting the metric field in this 
VJCF by (jfj^i, one has Lvac = f{R)\/—g- One makes the standard transformation to the Einstein frame, 
p = f'{R) and — pg^iv Solving the latter equation for R one finds R = r{p) and this defines the direct 
and inverse Legendre maps, p = and r = r{p). The field p, in the "phase space" of the system, becomes 
an independent dynamical variable and the relation p = f'{R) is a consequence of the equations of motion 
generated by the Helmholtz Lagrangian, or equivalently (after conformal rescaling) by the Einstein-frame 
Lagrangian (1.5), 

3 
V 

rip) f[rip)] 



with 



R 



2V{p) = 



(CIS) 
(C.14) 



p p 

Assuming that the rescaled metric g^j,,, is physical, one now adds the minimal-coupling term 2£jnat{g] ^') to 
the Lagrangian (C.13). This clearly affects the equations of motion and therefore the relation p = f'{R), 
holding in vacuum, fails to be true in the presence of matter. 

The present formalism restricts the allowed matter Lagrangians to those which do not depend on deriva- 
tives of the metric. Although stringent in itself, the restriction admits the physically important cases of 
perfect fluid, pure radiation, flelds of spin and usual fields of spin 1. The full Lagrangian is then 



R - ^r'P^f.P^u - 2V{p) + 2Cat(5; *) 



(ai5) 



The usual conformal rescaling ^^^1/ = P9iJ.v transforms (C.15) into the following Lagrangian for g^j^^, P and 
the matter: 

LAP-, 9,^) = V^{pR-2p''[V{p)-l^^,{pg;^)]} . (CIS) 
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It has the form of a Hc;hiiholtz Lagrangian (C.5) because th(^ transformation cancched the kinetic term for 
the scalar. The information about the original vacuum Lagrangian Lvac is encoded in the potential V(p). 
The requirement that does not depend on derivatives of p is met provided ^mat contains no covariant 
derivatives (thus, it may contain at most first derivatives of ^f). The Hamiltonian in the present case is 



H{p; g,^)=p- r{p) - f[r{p)] - 2pH^,t{pg; 



[C.n) 



and the equation for p becomes (2.16), R{g) — r{p) + g^'^Tj^v{^,P9) = ; the matter stress- energy tensor 
is defined, as always, in terms of the physical metric g^j^^ and expressed in terms of g^^. If ^,nat = or 
T^iy is traceless (e.g. matter is conformally invariant), then R{g) = r{p) and p — f'{R) will be a solution; 
otherwise the solution P{R;g,'^) ^ f'{R)- The Helmholtz Lagrangian (C.12) corresponds to a nonlinear 
Lagrangian (C.9) for g^^ and ^P. The NLG Lagrangian L^^^ so obtained reduces to the original vacuum one 
if the matter fields are set to vanish; however, in general it is not equal to the sum of Lvac and a matter 
term analogous to £ma.t (see the examples below) . In the purely metric picture, the relation between the two 
metrics is expressed by — P{R;g, 4')g^,y, which explicitly depends on matter. For this reason the MJCF 
metric corresponding to a given physical metric g^i, does not coincide in general with the VJCF metric ^j^jy, 
which is implicitly defined by g^v = f'{R)gij,v In MJCF the equations of motion for matter, (C.ll), read 



_5_ 
5* 



-gp^f^ma.t{pg; *) 







(C.18) 



The fourth-order field equations for gravity, (C.12), can be put in a number of equivalent forms. Here wc 
recast them in the form which is closest to that for the vacuum case. Upon using r{p) = R{g) + g^^Tfj^^ one 
finds that PGfj,v + \Pg^ivT{P) = PR^v + ^Pgij.ug"^Taf3, then defining a scalar M{R;g,^) as the matter 
contribution to P, P = f'{R) + M{R;g, one arrives at the following gravitational field equations: 



Q^^ + M{R; g, ^)R^^ - DM- - / Wl = 

= P{R; g, *)[V(P<?; *) - ^g^.g^^^T^^] 



(C.19) 



Here = /'(_R)i?^^ — ^f{R)g^i, — V fj.'V^f'{R) + .g^,y □ /'(i?) is the l.h.s. of equation (1.2) and satisfies 
the generalized Bianchi identity VQfj,v = (Appendix A). For traceless matter, (C.19) reduces to equation 
(2.13). One sees that matter and gravitational variables are inextricably intertwined in (C.19), thus these 
equations do not provide conservation laws similar to those in General Relativity. Formally, four matter 
conservation laws arise upon taking the divergence of (C.19), then the purely gravitational part, Q^i^, dis- 
appears and one is left with interaction terms. The resulting equations, however, are too complicated to be 
of any practical use even in the simples cases. 

Below wc give some examples of finding the nonlinear Lagrangian for a given form f{R) of the vacuum 
Lagrangian and for dust or a scalar field + electromagnetic field as the matter. 



Example 1 : quadratic Lagrangian and charged scalar field. 

Let the vacuum Lagrangian in VJCF be Lvac = {aR^ + R)\^—g ■ Then the "vacuum" inverse Lcgendre 
map is r(p) = ~ 1) ^^nd f[r{p)] = j^{p^ — '^)- In Einstein-frame metric g^^ = pg^v one adds the 
interaction with a massive complex-valued scalar field V minimally coupled to electromagnetic field. 



(C.20) 



Here D/j^ip = dutp — ieA^xp. After conformal rescaling ^^,y = p^^i/ this Lagrangian becomes (up to a full 
divergence) 



PR-^a^P- 1)' - P9^''D^xl,{D,xl,r - mVi>r - ^F^^.F^^g'^-'g^f 



'-9 



(C.21) 
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and the equation (C.14) is solved to yield 



P{R,g;ip) = 2 Am^iptp* + 



2a 



The nonlinear Lagrangian (C.9) in MJCF, 



(C.22) 



— F Fa o^''a°''^ — 



1 

2a 



generates the following fourth-order gravitational field equations: 

P{R, g; ^)R^^ - V^V.P + ff^. □ P - V-)* 



n 2 



2a 



4a 



(C.23) 



(C.24) 



Example 2: quadratic Lagrangian and dust. 

As previously wc choose in VJCF Lvac = {aR^ + R)^/—g , so the functions / [^(p)] and V{p) arc 

as in Example 1. In the Einstein frame one adds pressureless dust with energy density p, thus the full 
Lagrangian is 

3 



R 



ip-iY 



Upon conformal rescaling one finds the Helmholtz Lagrangian 



-9 



pR-j^{p-lf + 2p''p 



and accordingly the Legendre map in the presence of matter is 
Simple manipulations yield the nonlinear Lagrangian in MJCF, 



-^NL — 



4a 



(12a/5+l) 



2ai?+ 1 
2ap+ 1 



and the field equations for the metric, 



1 



iC.25) 



(C.26) 



(C.27) 



(C.28) 



P{R, p)R^, - ^i^VuP + 5/.. □ P - -^(P - l)gy.. = P'p{u^,u, + -g^^) 



(C.29) 



8a' /c^p^ . 2' 

It is very difficult to extract from these equations any nontrivial information about motion of dust particles. 
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Example 3: logarithmic Lagrangian and dust. 



Cases where the equation (C.18) for P can be explicitly solved are quite exceptional. Besides the 
quadratic Lagrangian and a simple matter source, one of these is provided by a logarithmic Lagrangian (an 

exponential one is not the case) and dust. 

In vacuum Jordan frame one has Lvac = ^ + <iR)\/~9 (for simplicity, we consider only the sector 
on which (1 + aR) > 0) and proceeding as in the previous cases one finds 



p = > , r{p) = -[--l) , f[r{p)] =--lnp 

ar + 1 a \p J a 



and 

The corresponding Legendre map is 



-(lnp+l)-l 
P 



pR- -{\np-p+l) + 2p'^p 



(CM) 



'r^'-)\'-^ 

, a J a 



a 



(there is another solution P{R) < which is discarded) and the resulting nonlinear Lagrangian is 




, a J a \ a 



a J a \ a 



+ 



-^[21n(8p)-l] 



(C.32) 



notice that L^^ = i ln(l + aR)sf^ + 2p(l + aRy^^f^ + 0{p^). 



Appendix D: The inverse problem of nonlinear gravity 

By means of the Helmholtz Lagrangian method and the conformal rescaling it is possible to map a 
nonlinear vacuum gravity theory with Lagrangian Lvac = \/^9 1{R) to a dynamically equivalent system 
consisting of Einstein gravity and a minimally coupled nonlinear scalar field with a potential determined 
by the function /. The inverse problem of nonlinear gravity consists in making an inverse transformation. 
Given a scalar field which self-interacts via an arbitrary potential U{(j)), is it possible to map the Lagrangian 



-9 R 



5^>,/x</',. - 2C/(0) 



(D.l) 



in the Einstein frame to an equivalent Lagrangian L^^ = f{R) in the VJCF? With one exception, the 
answer is "yes", but in most cases the nonlinear Lagrangian cannot be expressed int terms of elementary 
functions (even if U{^) can be). The procedure is as presented in Appendix C. First one redefines the scalar 
by setting <p = y/Jlnp and then makes the conformal rescaling g^^ = p^f^^ to cancel the kinetic term for the 
field p. After discarding a divergence term, {D.l) becomes a Helmholtz Lagrangian, 



L^=^\pR-2p^U{p)] , 



{D.2) 
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where U{p) = [/(yi'lnp). According to (C.7), the inverse Legendre map p = P{R) is a solution of the 
algebraic (w.r. to p) equation 

i?(ff)-^[2/C/(p)]=0; {D.3) 

then the function / is determined by the relation f'{R) = P{R)- 

An alternative method consists in solving an ordinary differential equation directly for /. Comparing 
(-D.2) with (C.16) one sees that U{p) = V{p), where V is given in (C.14) with p = ^^^^ and r{p) its inverse 
function. Then the equation U{p) = V{p) becomes a differential equation^^for /(r) 

■^{f'fUif) =rf'-f {DA) 

Solving this equation, however, is by no means easier than solving {D.3). Moreover, to ensure the consistency 

of the Legendre transformation, a solution to (£'.4) should also meet the R regularity condition ^ 0. 
Differentiating (-D.4) w.r. to r and using /" ^ one arrives at an equivalent equation, which added to (-D.4) 
yields 



dp 



thus / as a function of p is given by 



f[r{p)] = 2p^ 
while r{p) is determined from the form of V, 



U{P)+P^ 



{D.6) 



r{p) = 2pU{p) + i/[r(p)] . {D.7) 

To determine f{r) without any integration one should compute f[r{p)] from {D.6), then solve {D.7) for p{r) 
(which is the same as solving {D.3)) and insert the solution back to f[r{p)]. 

First we note that the method does not apply to the simplest case of the massless linear field, U{(j>) = 
0. In fact, in this case (-D.4) is solved by / = Cr and this solution is excluded because it is not ir- 
regular. Thus, the above- mentioned exception is provided by the Einstein-frame Lagrangian L{g, (f>) = 

\J —g R — g^'^(l).fj.<i',v , for which the (p-&eld cannot be eliminated to yield an equivalent purely metric NLG 
theory. This Lagrangian, on the other hand, can be turned into a STG Lagrangian (2.18) by a suitable 
conformal rescaling. A constant potential in {D.l) is interpreted as a cosmological constant, (/{(p) = A. 
Then {D.6) and {D.7) easily yield f{R) = -^R^- This result is also obtained as a particular ("singular") 
solution of (J3-4), which in this case is a Clairaut equation (and has also a general solution, which should be 
excluded as being linear). 

For nonconstant potentials U{4>) the solutions do exist, but these arc practically inaccessible since one 
is unable to solve {D.3) or {D.7). For the most interesting from the field theoretical viewpoint - potential, 
U{cj)) = A(/)", A = const., n = 2,3,4, one finds U{p) = ^i{\np)'\ with /.i = (|)^ A, 

and f[r{p)] = 2^^p^{\np + n){lnp)"-^ {D.8) 

r{p) ^ 2^ip{2 \tlp + n)(lnp)"-i . (£).9) 

The Liouville field theory [32] provides one of the few examples where equation {D.7) can be explicitly 
solved. In this case U{4>) = Ae"'^', {A, a constants), and one finds f[r{p)] = 2(1 + /3)Ap'^+^, with P = -/fa; 
r{p) = 2(2 + l3)Ap^+'^, then 

and P{R) = [2{2 + I3)A]~'^R^ 

+ 2 f3 + 2 

f{R) = 2(1 + p)A[2{2 + P)A] ^ i?^ . {D.W) 



In [52], Teyssandier and Tourrenc present a system of differential equations which is equivalent to {D.A). 
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Appendix E: Local dynamical equivalence with GR and global solutions of NLG theories 

In Appendix C we have reviewed the Legendre-transformation method. As in classical mechanics, the 

relationship between the original fourth-order dynamics and the second order equations generated by the 
Helmholtz Lagrangian is a "dynamical equivalence" . This means that there is a one-to-one correspondence 
between solutions of the two systems of equations. This correspondence is provided by the Legendre map. 
Since the Legendre map is not. in general, globally invertible, the correspondence (and thus the dynamical 
equivalence) holds locally. A distinct problem is whether the conformal rescaling leading from the Jordan 
frame to the Einstein frame can be defined globally or not. In the case of quadratic Lagrangians only 
the second problem arises. If, for a given NLG lagrangian (1.1), both the inverse of the Legendre map 
and the conformal rescaling can be globally defined for all solutions of the NLG equations (1.2), then the 
set of these solutions actually coincides with the set of solutions of the Einstein-frame equations (General 
Relativity + scalar field). Both conditions, however, do not hold in general, and therefore one should make 
a distinction between local solutions and global solutions. This does not mean that the local dynamical 
equivalence with GR is irrelevant while dealing with global solutions of NLG; on the contrary, it allows to 
understand completely how global solutions of NLG theory look like. In this Appendix we discuss separately 
the problems arising when the Legendre map is not bijective and when the conformal factor vanishes. 

Once one is given a NLG lagrangian (1.1), the Legendre map is defined, as it is explained in Appendix 
G, by the function p = f'{R). Whenever the Legendre map can be inverted, one constructs the Helmholtz 
Lagrangian and the corresponding second-order equations. The dynamical equivalence of the latter with 
NLG equations means, in one direction, that for any solution of the second-order system (C.3), represented 
in a local coordinate system by {gij,,^{x°'),p{x°')), the metric tensor g^ivix"^) is a solution of the original fourth- 
order NLG equations. Conversely, to any given solution .9^1^(2;") of the NLG equations, one can associate the 
corresponding curvature scalar R = i?(.T"); one introduces the scalar field p(x") — /'[i?(x")], and the pair 
(gfj,i^{x"),p{x")) then fulfills the system (C.3). The notion of dynamical equivalence means nothing more 
than that. In particular, the value of the Helmholtz Lagrangian for a given pair {g,p) coincides with the 
value of the NLG Lagrangian (1.1) for the same metric g only if {g,p) is a solution of the field equations; two 
dynamically equivalent action functionals have the same stationary points, but they take different values at 
states of the system which do not correspond to classical solutions. 

The dynamical equivalence holds under the i?-regularity condition f"{R) ^ 0. For the sake of simplicity, 
we first consider a vacuum NLG (1.1) with / analytic. Then the equation 

/"(i?)=0 {E.l) 

is an algebraic equation which has a discrete set of solutions R — ki (the label i runs over a finite set of 
integers if / is polynomial). To these values of the curvature scalar corresponds, under the Legendre map, 
a set of values p = Ci = p{ki) of the scalar field (the correspondence between the sets {hi} and {q} is not 
bijective, in general). The equations of the system (C.3) are not defined for p — Ci, because the domain 
of the inverse Legendre map r(p), which occurs explicitly in (C.3), does not include these points. A global 
dynamical equivalence holds whenever one can ensure that, for all possible solutions of the NLG equations, 
the curvature scalar nowhere attains one of the values ki at any point of space time. This is actually possible 
only if {D.l) has no real solutions at all: this happens, e.g., for J{R) — R + aR?, e"^, cosh(ai?), OkR^^ 
with a/j > 0. 

If the set of critical curvatures {ki} is not empty, the solutions of NLG equations (1.2) can be divided 
into three groups. First, the NLG equations may admit constant curvature solutions for which R{x°') = ki 
(for some i) in the whole space-time. Such solutions do not correspond to any solution of the second-order 
system (C.3) . However, in both cases it is very easy to check by direct computation whether such solutions 
exist, and take them properly into account. The second group of solutions includes global solutions which 
do not attain any critical value at any point; for each solution of this type we can find a corresponding 
global solution of (C.3). The third possibility is that a (non-constant) solution of NLG equations attains 
some of the critical values in some regions of space-time. Under reasonable assumptions on the regularity 
of the metric tensor, this region has measure zero: it is a lower-dimensional submanifold of space-time. 
This critical submanifold, which is different for each solution, may have several disconnected components, 
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corresponding to different critical valnes. The critical submanifold separates space time into open domains 
on which the i?-regularity condition is fulfilled. Inside each regularity domain, one can associate to the global 
NLG solution under consideration a local solution of the system (C.3). Such a local solution (51^1,(0;"), 
has the property that p{x") Ci, for some i, when one approaches the critical submanifold, thus it can 
be extended by continuity across the critical surfaces. As a matter of fact, the system (C.3) itself can be 
extended to such points as well: the derivative of r{p) blows up at critical points, but r{p) admits a finite 
limit there. 

The most relevant problem which can arise from a failure of the regularity condition on lower dimensional 
submanifolds is connected with the uniqueness properties of solutions. A given local solution may admit 
several continuations outside its regularity domain. This is connected with (but not entirely dependent on) 
the following fact, which we have not yet taken into account: whenever f'{R) is not bijective, then in each 

regularity domain the inverse function may not be unique. 

In general, the equation f'{R)\j^_^^^^ — p = admits several roots for a given value of p. Assume for 
instance that f'{r) is a polynomial of third degree; then, for a generic choice of the coefficients there are 
two critical values {ki,k2} (fci < ^2), and one finds one root ri{p) for —00 < p < ki, three distinct roots 
{^i{p)}i=2,3,4 for k\ < p < k^ and again one root r^{p) for p > k2- Putting together the families of roots 
which have the same limit for p ^ ki, one obtains exactly three inverse Legendre maps: r_(p), defined for 
p G (—00, ^2); ^0(^)7 defined only for p € (fci, fc2); and r^{p) for p G (fci, 00). Thus one has to introduce three 
distinct Helmholtz Lagrangians, each one containing a different potential for p. Each Helmholtz Lagrangian 
is defined only for the range of values of p for which the potential is defined. In the example considered 
above, for p G (fci, fc2) there are three different second-order scalar tensor systems (C.3) locally equivalent to 
(1.2) (we may then speak of different "sectors" for the scalar field). It might seem that this causes troubles 
such as ill- posedness of the Cauchy problem, but it is not so. In fact, the ranges for p may overlap (and 
they do, in general), but the ranges of the curvature scalar R corresponding to different potential are always 
disjoint, since the Legendre map p = p{R) is globally and uniquely defined. Thus, if one considers a global 
NLG solution, in each regularity domain one knows exactly to which sector the scalar field belongs, just by 
checking the range of values of i?, and therefore one knows without ambiguity which scalar tensor equations 
(C.3) are fulfilled by the pair {g,p = f'{Rj). If, in the opposite direction, one considers a Cauchy problem 
for the scalar-tensor model, then the Cauchy data are compatible with at most one sector of the theory, and 
therefore determine the equations which have to be solved. 

Hence, the local dynamical equivalence which holds in the cases we arc dealing with can be described as 
follows: to each solution of the NLG equations (1.2), with the exception of the solutions for which f"{R) = 0, 
corresponds a scalar-tensor pair {g,p). This pair has one or more i?-regularity domains, which are the open 
regions of space time on which f"{R) 7^ 0. Inside each regularity domain, the pair {g,p) coincides with 
a solution of the scalar-tensor model described by a suitable Helmholtz Lagrangian (C.2), where r{p) is 
the unique inverse Legendre map defined for the range of values of R occurring in that regularity domain. 
Conversely, all solutions of NLG equations (besides the particular constant-curvature solutions mentioned 
above) can be recovered by matching together local solutions of the scalar-tensor models corresponding to all 
possible inverse Legendre maps r{p) (for the given NLG Lagrangian); each local solution exists, by definition, 
only in the open region of space time on which the scalar-tensor equation (C.3) is well-defined; however, in 
general the solution can be extended to the closure of the domain of definition of the equations, and can 
therefore be matched with analogous local solutions on adjacent domains; the minimal glueing prescription is 
that the curvature scalar of the metric tend to the same limit on both sides of a critical surface. This entails 
the continuity of p, but is a stronger assumption because R could have a discontinuity, due to the jump 
from one sector to another, even when p is continuous across the critical surface. The matching requirement, 
therefore, prevents certain pairs of sectors to occur in adjacent domains. 

The existence of different sectors depends on the non-uniqueness of the inverse of f'{R)- By a different 
parametrization of the scalar-tensor model, namely by using the alternative Helmholtz Lagrangian (C.2'), one 
avoids this problem. The Helmholtz Lagrangian, and therefore the scalar-tensor equations, need nevertheless 
to be continued to critical surfaces (whereby the equation (1.2) cannot be recast in normal form); the 
continuation may be not regular enough to avoid well-posedness problems, and the distinction between 
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global and local dynamical equivalence should be kept. Moreover, we have already mentioned in Appendix 
C that the definition of the physical scalar field requires in any case the use of an inverse Legendre map; 
thus, the physical picture of the theory does include different sectors, unless the NLG Lagrangian is globally 
i?-regular. 

In Appendix C, we have also showed how a nonlinear Lagrangian can be recovered, by inverse Legendre 
transformation, after a modification of the Helmholtz Lagrangian by the addition of matter coupling terms; 
however, we did not take into account the possible non-uniqueness of the Helmholtz Lagrangian. We now 
complete the disciission on this point by checking that, in the case of multiple Helmholtz Lagrangians 
corresponding to the same vacuum NLG Lagrangian, matter coupling terms can be consistently added in all 
sectors, so that the inverse Legendre transformation leads to a unique nonlinear interaction Lagrangian. More 
explicitly, suppose that for a given f{R) in the NLG Lagrangian(l.l) there exist several inverse Legendre 
maps, {ri{p)}, defined on different ranges of the scalar field p (which partially overlap). Accordingly, one 
has several locally equivalent Helmholtz Lagrangians, each one containing a different "Hamiltonian" (we 
follow the notation of Appendix C) Hi{p) = p ■ ri{p) — f[ri(p)]. Now, what happens if one modifies each 
Helmholtz Lagrangian as described by equation (C.16)? One would reasonably expect that matter couples 
to the metric in the same way for all sectors, since in the Einstein frame matter is not coupled at all with 
the scalar field p: the "Hamiltonians" should hence become Hi{p) = p ■ ri{p) — f\ri{p)] — 2p^i'mat(p.9; V')- Do 
such "Hamiltonians" correspond to a unique nonlinear Lagrangian? Using the explicit form of Hi[p), the 
equation (C.7) becomes 

R{9) - nip) + 4pCat(w; V-) + 2p^j^9^u = {E.2) 

Let us redefine the auxiliary scalar field by setting p = f'{u), as explained in Appendix C: one has in each 
sector u = ri{p), i.e. u is a multivalued function of p; for u, however, the following unique equation holds: 

R{g) -u + 4/'(n) Cat[/'(«)5; V'] + 2[f'{u)]^^g^. = (£.3) 

leading to solutions u — U{R;g,^) which are independent of i. For any such solution, P{R;g,ip) = 
f'[U{R;g,4>)] solves (-D.2) for all i, wherever the domain of ri{p) intersects the range of P{R;g,tjj); thus, 
we see that even for ranges of values of p for which several Helmholtz Lagrangians coexist, {E.2) admits 
solutions which do not depend on i. A nonlinear Lagrangian can then be defined by setting 

L,, = P{R;g,nR{9)-U{R;gA)]V^ + f{U)V^ +2[fiU)fe^MU)9-,i^]V^ 

For a more detailed investigation of this topic, we refer the reader to [53]. 

We now consider the second problem: namely, whether the Einstein conformal frame corresponding to 
a solution of NLG equations can be globally defined or not. If we assume, as we did in Section 1, that the 
conformal factor coincides with the scalar field p = f'{R), the Einstein frame is globally defined provided 

f'{R) > {EA) 

everywhere. Otherwise, the rescaling introduces new singularities. This problem has been raised by several 
authors; most of them, however, consider the procedure leading from the metric g in (1.1) to the pair (g, <p) 
occurring in the Einstein-frame Lagrangian (1.6) as a single transformation, while it actually consists of three 
distinct steps (addition of a new independent variable, conformal rescaling, redefinition of the scalar field 
(f) ~ ln(p)). For instance, Maeda in [5] comments about the local character of the equivalence between the 
two frames, but he refers only to possible failures of condition (i?.4), although the proof of the equivalence 
which he gives in the same paper is valid only if {E.l) is fulfilled; moreover, his comment refers mainly to 
the problem of defining the scalar field (j) (denoted by ij) in his paper) at points where f'{R) = 0, rather than 
to the fact that the Einstein metric vanishes at these points. 

The transformation leading from the scalar-tensor equations (C.3) to the Einstein-frame equations (1.4), 
which is the only genuine conformal rescaling occurring in this context, is mathematically a mere change of 
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variables. In contrast to the case; of the; dynamical cqnivalcncc relating (1.1) to the Helmholtz Lagrangian, 
the latter is now simply re-expressed in terms of the new variables, so both the Helmholtz Lagrangian and 
the Einstein-frame Lagrangian (1.5) take the same values (more precisely, differ by a total divergence) at 
conformally-related pairs {g,p) and {g,p), even if these pairs arc not solutions of the field equations (this 
fact is mentioned by Wands [6] for the case of conformal equivalence between JBD theory and GR+scalar 
field). The condition {EA) is, in principle, completely independent from {E.l). Only for particular NLG 
Lagrangians, e.g. ,f{R) = R^, with A; > 2, p = is also a critical value for the Legendre map, and the two 
conditions partially overlap. 

In some cases, we can ensure from the very beginning that p is positive everywhere. If f{R) is a 
polynomial of odd degree with positive highest-order coefficient, then the Legendre map p{R) is bounded 
fom below, and the coefficients of the lower-order terms in f{R) can be chosen so that p > Pmin > 0. 

If f{R) is a polynomial of even degree (for instance, in the quadratic case), then the scalar field p will 
range over the full real line. In this case, we should again regard the solutions of NLG as divided into three 
groups. For some solutions we may have p = 0; such solutions have constant scalar curvature and can be 
easily singled out and treated separately (under the assumption (3.8), flat space dose not belong to this 
class of solutions). For other solutions, {EA) holds everywhere in space-time and the Einstein frame is then 
globally defined. The problem arises for the third group of solutions, for which p vanishes on some region of 
space-time. As previously, under suitable assumptions such regions have zero measure and in most cases are 
lower-dimensional submanifolds of space-time. For these solutions, the conformal transformation introduces 
singularity surfaces. In rigorous mathematical terms we should say that the equivalence holds only where 
both and g^j^ are nondegenerate; both the Lagrangian and the field equations for g^i, are defined only 
where g^^ is regular. However, for some Lagrangians the Einstein-frame equations admit a solution (g,p) 
such that Qfj,,^ is nondegenerate but p vanishes, and therefore it is the Jordan-frame metric which cannot be 
recovered at such points: whether we should regard this situation as "singular" or not, from the physical 
viewpoint, depends on the physical significance attached to the metric. 

A separate problem, which has an evident physical relevance, concerns the signature of the rescaled 
metric. Whenever p < 0, the conformal rescaling can be defined but it changes the signature of the metric. 
If we assume that the Einstein-frame metric is the physical one, this is not acceptable because the choice 
of the signature, although arbitrary, should be globally consistent. This problem is solved by some authors 
(e.g. Barrow & Cotsakis and Maeda [5]) by changing the definition of the Einstein metric to 



This introduces a discontinuity in the first derivatives of the metric, but this happens exactly on the singu- 
larity surfaces, and constitutes a minor problem. The definition, however, causes an overall minus sign to 
appear in front of the Einstein-frame Lagrangian when p is negative (in other words, the; entire Lagrangian 
is multiplied by -pr; this seems to make the Lagrangian discontinuous at p = 0, but in fact the Lagrangian 



An alternative viewpoint is to reject the requirement that both metrics have the same signature. Since 
we claim that only one metric is physical, any requirement on the signature should affect only this one. Once 
we have found a physically acceptable solution {g,p), the possible fact that the corresponding Jordan-frame 
metric changes signature due to a change of sign of p may be regarded as irrelevant (anyhow, to change 
signature, the Jordan frame metric should necessarily become singular somewere, and this feature is much 
more relevant). 



9tiu = \p\gt,,.- 



{E.5) 
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Fig. 1: Relations among NLG theories (the Lagrangian (1.1) is expanded in power series around R = 0), 
STG theories and General Relativity. The connection (A) is described in Appendix C, (C.1,2); (B) is 
represented by (2.21,22); (C) is outlined in Sect. 1 (1.3,6); (D) and (E) are discussed at the end of Sect. 2; 
(F) represents mappings between STG theories with different co{(p) (and different A(<^), if any), being a 
combination of an arbitrary conformal rescaling and a suitable field redefinition. 
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